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Abstract 

The aim of this paper is to rigorously investigate the orbital magnetism of core electrons 
in crystalline ordered solids and in the zero-temperature regime. To achieve that, we consider 
a non-interacting Fermi gas subjected to an external periodic potential within the framework 
of the tight-binding approximation (i.e. when the distance R between two consecutive ions is 
large). For a fixed number of particles in the Wigner-Seitz cell and in the zero-temperature 
limit, we write down an asymptotic expansion for the bulk zero- field orbital susceptibility 
and prove that the leading term is the superposition of the Larmor diamagnetic contribution 
(reducing to the Langevin formula in the classical limit) together with the so-called orbital 
Van Vleck contribution. 



Contents 

1 Introduction & the main result. 

1.1 An historical review-Introduction. 

1.2 The setting and the main result. 

1.3 Discussion on the assumptions. . 

1.4 An open problem 

1.5 The content of the paper 



Appendix. 

4.1 Proof of Propositions 13.61 and 13.71 . 

4.2 Proof of Propositions [SH and [STOl 

4.3 Proof of Lemma 13.131 

4.4 Proof of Lemmas li?34.5l 



2 An approximation of the resolvent via a geometric perturbation theory. 

3 Proof of Theorem [HB |l2 

3.1 Proof of (i) n2 

3.1.1 The location of the Fermi energy 12 

3.1.2 Isolating the main contribution at zero-temperature 15 

3.1.3 Isolating the main contribution at zero-temperature - Continuation and end. 17 

3.2 Proof of (h) 18 

3.3 Proof of (iii) [22 



5 Acknowledgments. l33l 

PACS-2010 number: 75.20.-g, 51.60. +a, 75.20.Ck, 75.30.Cr 
MSC-2010 number: 81Q10, 82B10, 82B21, 82D05, 82D20, 82D40 

Keywords: Orbital magnetism; Atomic magnetism; Zero-field susceptibility; Langevin formula; 
Larmor diamagnetism; Van Vleck paramagnetism; Tight-binding approximation; Geometric per- 
turbation theory; Gauge invariant magnetic perturbation theory. 



*Department of Mathematical Sciences, University of Aarhus, Ny Munkegade, Building 1530, DK-8000 Aarhus 
C, Denmark; e-mail: baptiste.savoie@gmail.com 



1 Introduction &; the main result. 



1.1 An historical review— Introduction. 

The first important contribution to the understanding of diamagnetisni of ions (we assimilate 
an atom to an ion of charge zero) and molecules (including the polar ones) go back at least to 
1905 with the three papers [5^1133133 of P. Langevin. We mention all the same that, in all likeli- 
hood, W. Weber brought between 1852 and 1871 the pioneer ideas through his molecular theory 
of magnetism in which he already introduced the idea that the magnetic effects are due to orbiting 
motion of electric charges around fixed charges of opposite sign, see [4]. The Langevin's theory 
essentially leans on the classical Maxwell equations of electromagnetism. Putting things back into 
context (the atomic structure was experimentally discovered in 1911), Langevin considered that 
matter is formed by electrons in stable periodic motion (the mechanical stability being ensured 
by the mutual actions between electrons); in particular, the molecules contain at least one closed 
electron orbit with a fixed magnetic moment out of the field (electrons are assimilated with partic- 
ulate Ampere's currents), and the different orbits in each molecule have such a moment and such 
orientations that their resultant moment may vanish, or not. Starting with these assumptions, 
he calculated the mean variation of the magnetic moment (orthogonal to the orbit) of electron 
moving in intramolecular closed orbits under the influence of an external constant magnetic field. 
This led to the so-called Langevin formula for the diamagnetic susceptibility per unit volume of 
n (intramolecular) electrons , see [38j pp. 89]: 

XLi-^n-^{r'), (1.1) 

where r is the distance from the molecule's centre of mass to the electron, considering its motion in 
the projected orbit on the plane perpendicular to the magnetic field; (r^) is the average over all 
the molecular electrons. The obtained result is independent of the temperature and independent 
of whether or not the initial resultant magnetic moment of the molecule is null. For the sake of 
completeness, we mention that Langevin derived also in |38j the analytic expression of the Curie's 
empirical law for molecules of paramagnetic substances: the paramagnetic susceptibility per unit 
volume of N identical molecules with a non zero resultant magnetic moment M reads as, see [351 
pp. 119]: 

xr" = fw (1.2) 

In 1920; W. Pauli was interested in diamagnetism of monoatomic gases in [32] ■ Still within 
the framework of the classical theory, his approach slightly differs from the one of Langevin in 
the sense that it is based on the Larmor's theorem. This latter states that the introduction of an 
external constant magnetic field with intensity B causes the angular velocity of an electron in a 
periodic orbit to be increased by ~Be/2mc without the orbit undergoing any modification. This 
led to the formula for the diamagnetic susceptibility per unit volume of N identical atom supposed 
to have random orientations in the space, see [HI pp. 203]: 

xs^^-^^E^> (1-3) 

i 

where the sum is over all the electrons of a single atom, is the square distance from the nuclei 
to the i-th electron and has to be understood as the time average of rf. The formula (|1.3p is 
sometimes referred as the Langevin formula in the form given by Pauli. 

In 1927-1928, J.H. Van Vleck revisited in a series of three papers [48l|49l[50] the Langevin the- 
ory on magnetic susceptibility (as well as the Debye theory on dielectric constant) of ions/molecules 
but within the framework of the 'new' quantum mechanics. Starting from the assumption that 
each ion/molecule has a permanent dipole moment of constant magnitude with slow precessions 
(i.e. the precession frequencies are small compared to 1//3), Van Vleck derived a general for- 
mula in [48l Eq. (13)] for the total magnetic susceptibility of N (randomly oriented) identical 
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ions/molecules in the zero-field limit. Note that only the contribution coming from the electrons 
is taken into account, the one coming from nuclei is considered as negligible. The formula is a 
generalization of the complete classical Langevin formula (obtained by adding (jl.3|) and (|1.2[) ) 
but including quantum effects, and consists of the sum of two contributions. The first one (it 
corresponds to the second term in the r.h.s. of |48l Eq. (13)]) is /3-dcpendent (the dependence 
is in 1//3), always paramagnetic and arises from the presence of a non-zero permanent magnetic 
moment; ergo it is identically zero when considering atoms since as a rule they have no permanent 
dipole moment. Moreover its classical equivalent is (|1.2|) . The second one is A^a, where a is a /?- 
independent constant representing the induced magnetic moment. It consists of the superposition 
of two terms, see HH Eq. (15)]: 

Na = 7V( JTvv + ^La). (1.4) 

Here ^i^a, is the so-called Larmor contribution which is purely diamagnetic, and N ^i^^i reduces 
to the diamagnetic Langevin formula in (jl.3p in the classical limit. As for St^wf., it is purely 
paramagnetic and has no equivalent in the classical theory (it will bear the name of Yan Vleck 
susceptibility later on). Furthermore, Van Vleck analyzed the origin of each one of these two 
contributions, see e.g. [SU Sec. VIII.49]. Restricting to the case of a single ion, he showed that 
^■vy and are generated respectively by the the linear part and the quadratic part of the 
Zeeman Hamiltonian of the electron gas which are defined by: 

-ffz, linear := A*s(L + 5oS) • B, Fz, quadra (^""i X B^ . (1.5) 

i 

Here L and S stands for the total electronic orbital and spin angular momentum respectively, go 
and fj,B for the electronic (7-factor and Bohr magneton respectively, B for the external constant 
magnetic field and for the position vector of the i-th electron. Formulated with the notations 
of the modern quantum mechanics, the contributions ^yy and in (|1.4|) are given in ()1.6|) . 
From the foregoing. Van Vleck concluded that ,^?^a always exists and is in competition with the 
paramagnetic contribution ^yy when the ion has either its total electronic spin angular momen- 
tum or orbital angular momentum different from zero in its 'normal state' (i.e. ground state). 
Moreover he claimed that ^yy does not vanish in great generality in the case of molecules, see 
[511 Sec. X.69]. Hence iV<^a is an upper bound limit to the diamagnetism of electrons in any 
case (ions/molecules). 

In regards to the method used by Van Vleck to derive the generalized complete Langevin formula, 
the starting point in [50l consists in using the Maxwell-Boltzmann statistics to express the 
magnetic susceptibility from the matrix elements of the average value of the magnetic moment 
per ion/molecule induced by an external constant magnetic field with intensity B; these matrix 
elements being defined as the first derivative w.r.t. i? of the energy in field. Next he used the 
stationary perturbation theory to expand in first powers in B (at least two) the energy in field, and 
then obtained such an expansion for the matrix elements of the magnetic moment. After expand- 
ing the exponentials Boltzmann distribution factor in power series of i?, the zero-field magnetic 
susceptibility is obtained by identifying the term proportional to B within the linear approxima- 
tion theory. The rest requires the permanence of moment and the slowness of precession postulates 
in order to isolate the /3-independent part from the /3-dependent one. 

Still within the framework of the quantum mechanics, another approach is commonly encoun- 
tered in Physics literature to treat the magnetism of ions/molecules, and does not involve (at 
least directly) statistical mechanics, see e.g. [Sj Chap. 31]. Provided that the intensity B of the 
external magnetic field is sufficiently weak, the method consists in computing the changes in the 
energy level of an ion by treating the magnetic field as a perturbation of the ground state energy. 
Expanding up to the second order in B the changes in energy level by the stationary perturbation 
theory, then the magnetic susceptibility is obtained by performing the second derivative w.r.t. B. 
This procedure can be justified from a physical viewpoint since the free energy is equal to the 
ground state at zero-temperature. Denoting by |0) the ground state and \n) n > the excited 
states of the ion, and assuming that B = Be^, k G {1,2,3} then the /3- independent part of the 
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magnetic susceptibility per unit volume of an ion is given by 



|(0|/4(L + .goS)-efc|n)p 




5](0|(r,xe..f|0). (1.6) 



We emphasize that the above formulae hold only in the absence of degeneracies. As pointed out 
by Van Vleck, only in the case when the ion has its total electronic spin angular momentum and 
orbital angular momentum both equal to zero, then the a reduces to the Larmor contribution 
^La- As a consequence of the Hund's rules, this is the case of an ion having all its electron shells 



The aim of this paper is to rigorously revisit the atomic orbital magnetism (by orbital mag- 
netism we mean that we focus only on the magnetic effects which not arise from spin effects). Our 
approach is substantially different from the ones mentioned above. To model the core electrons 
of an ion, wc consider a non-interacting Fermi gas subjected to an external periodic potential 
within the tight-binding approximation. Under this approximation, we suppose that the distance 
between two consecutive ions is sufficiently large so that the Fermi gas 'feels' mainly the potential 
energy generated by one single nucleus. For a fixed number of particles in the Wigner-Seitz cell, 
we write down an asymptotic expansion for the bulk zero-field orbital susceptibility of the Fermi 
gas under the tight-binding approximation and in the zero-temperature limit. The leading term 
of this expansion plays the role of the a. We emphasize that the bulk zero-field susceptibility is 
derived from the usual rules of the quantum statistical mechanics. The rest consists in isolating 
and identifying the Larmor contribution from the Van Vleck contribution. Unlike the above men- 
tioned works, it is found that the Van Vleck contribution has a positive sign when the number of 
particle is equal to one, otherwise it can be written as a sum of a positive and negative term. 

In the frame of the mathematical-physics, the rigorous study of orbital magnetism and more 
generally of diamagnetism, was the subject of numerous works. Let us list the main ones among 
them. The first rigorous proof of the Landau susceptibility formula for free electron gases came as 
late as 1975, due to Angelescu et al. in [1]. Then in 1990, Helffer et al. developed for the first time 
in [3D] a rigorous theory based on the Peierls substitution and considered the connection with the 
de Haas- Van Alphen effect. These and many more results were reviewed in 1991 by Nenciu in [40] . 
In 2001, Combescure et al. recovered the Landau susceptibility formula in the semiclassical limit 
in [13]. In 2012, Briet et al. gave for the first time a rigorous justification of the Landau-Peierls 
approximation for the susceptibility of Bloch electron gases. Finally we mention the following 
papers [231 1241 125j in connection with atomic magnetism. 

1.2 The setting and the main result. 

Consider a 3-dimensional quantum gas composed of a large number of non-relativistic identical 
particles, with charge (7 7^ and mass m = 1, obeying the Fermi-Dirac statistics, and subjected to 
an external constant magnetic field. The particles possess an orbital and spin magnetic moment. 
Since we only are interested in orbital effects and we do not take into account the spin-orbit 
coupling, then we disregard the spin of particles. Moreover, each particle interacts with an external 
periodic electric potential modeling the ideal lattice of fixed ions in crystalline ordered solids. 
Furthermore the interactions between particles are neglected (strongly diluted gas assumption) 
and the gas is at equilibrium with a thermal and particles bath. 

Let us make our assumptions more precise. The gas is confined in a cubic box centered at the 
origin of coordinates given by Kl := (— L/2, L > 0; we denote by |Al| its Lebesgue-measure. 

We consider a uniform magnetic field B := (0, 0, B), parallel to the third direction of the canonical 
basis of R'^, and we use the symmetric (transverse) gauge, i.e. the magnetic vector potential is 
defined by i?a(x) := x x = (— 2:2, xi, 0), 63 (0,0,1). In the following we denote by 
h := qB/c e M the cyclotron frequency. The potential energy modeling the interaction between 
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each particle and the ideal lattice of fixed ions is given by: 



Vr <■ - ^^)' ^ > 0, (f .7) 

where the single-site potential u satisfies the following assumption: 
{s^-c) u e C"'^(R'^;R) is compactly supported. 

We denote by fifl; the Wigncr-Scitz cell of the i?Z'^-lattice centered at the origin of coordinates. 
Introduce now the one-particle Haniiltonian. On C^(Ai) define Vi? > the family of operators: 

HR^L{b):^\{-iV -bsif + Vr, 6eM. (1.8) 

It is well-known that Vi? > and Vfo G M, (|1.8|) extends to a family of self-adjoint and semi-bounded 
operators for any L e (0, oo), denoted again by HR^L{b), with domain 'HJ(Al) n 'H^(Al). This 
definition corresponds to choose Dirichlet boundary conditions on dKi^. Moreover, by standard 
arguments Hrj^IJ}) has compact resolvent, and its spectrum is purely discrete with an accumu- 
lation point at infinity. We denote by {A^^^(6)}j>i the set of eigenvalues of i^i^,L(fe) coimting 
multiplicities and in increasing order. As well, Vi? > and V6 € K denote by Nr^l{E), E eM. the 
number of eigenvalues (counting multiplicities) of the operator HR^L{b) in the interval (— oo,i?). 
When Al fills the whole space, on C^(]R'^) define Vi? > the family of operators: 

HR{b):^^{-iV -bsi)^ +Vr, beR. (1.9) 

By [131 Thm. X.22], then Vi? > and V6 G K pT^ is essentially self-adjoint and its self-adjoint 
extension, denoted again by HR{b), is bounded from below. Sometimes we will use the shortcut 
notation Hr ~ HR{b = 0). Moreover, the operator HR{b) only has essential spectrum since it 
commutes with the usual magnetic translations. Further, Vi? > and V6 G M the integrated 
density of states of the operator HR{b) exists, see e.g. [311 Thm. 3.1], and it is given by the limit: 

Nr[E) :^ lim ^^f (f^^ ^ hm TrL^(«-)faA^P(-°°^f )(^^^^^ W)XA J ^ ^ ^ (l.^O) 
itoo |Al| itoo |Al| 

where XA^ denotes the multiplication operator by the characteristic function of A^, and Pi{HR^L{b)) 
the spectral projection of HR^Lib) corresponding to the interval i C M. 

Let us now define some quantities related to the Fermi gas introduced above within the frame- 
work of the quantum statistical mechanics. In the grand-canonical ensemble, let (^, z, |Al|) be the 
fixed external parameters. Here /? := (fcsT)^^ > {ks stands for the Boltzmann constant) is the 
'inverse' temperature and z := c^^ > (/i G K stands for the chemical potential) is the fugacity. 
For any /? > 0, z > and & G M, the finite-volume pressure and density are respectively defined 
Vi? > by, see e.g. [311 111]: 

1 -. oo 

PrAP, ^, b) := ^Tri.(A,){ln(l + ze-^^^-^))} = ^ ^ ln(l + ze-^'''^-^^% (1.11) 



L\ 



. (1,12) 

Let us note that the series in ()l.lip - ()1.12p are absolutely convergent since V6 G M and Vi? > the 
semigroup {e~'^^"-^(^), /? > 0} is trace-class, see e.g. [101 Eq. (2.12)]. Moreover, from [TOl Thm. 
1.1], then V/^ > and Vi? > Pr.l{I3, • , • ) is jointly real analytic in (z, 6) G x M. This allows 
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us to define the finite-volume orbital susceptibility as the second derivative of the pressure w.r.t. 
the intensity B of the magnetic field, see e.g. |2] and [H Prop. 2]: 



A'h,l(/3,z,6) ( ^ ) -^^{I3,z,b), /? > 0, z > 0, 6 G M, i? > 0. 



When Al fills the whole space, then the thermodynamic limits of the three grand-canonical quan- 
tities defined above generically exist, see e.g. [T2j Thms. 1.1 & 1.2] and [Til Sec. 3.1]. Denoting 
V/? > 0, Vz > 0, V6 e R and Vi? > the bulk pressure by Pr{/3, z, b) := limi^oo -Pr,l(/3, z, b), then 
we have the following pointwise convergences: 

Pfl(/3,z,6):=/3z^(/3,z,6)= lim/3z^^(/3,z,6), (1.13) 

and the limit commutes with the first derivative (resp. with the second derivative) of the pressure 
w.r.t. the fugacity z (resp. w.r.t. the cyclotron frequency 6). Moreover, V/3 > and Vi? > 
Pr(/3, • , • ) is jointly smooth in (z, 6) £ R!^ x R, see e.g. [i5] . 

Now assume that the density of particles pQ > becomes, in addition with the inverse temper- 
ature, a fixed external parameter (canonical conditions). Seeing the bulk density as a function of 
the ^-variable, denote Vii > by /i^''(/3, po, 6) G R the unique solution of the equation: 

PR{P,ti,b) ^ po, f3>0,beR,R>0. 

The inversion of the relation between the bulk density and the chemical potential is ensured by 
the fact that V/3 > 0, V6 € R and Vi? > 0, i—> /3_r(/3, /i, b) is a strictly increasing function on R, 
and actually it defines a C°°-diffeomorphism of R into (0, oo), see e.g. [151 [TT]. In the following 
we consider the situation in which the density of particles is given by: 

Po(i?)-^, noeN*,i?>0. (1.15) 
|"fi| 

Let us note that in (|1.15p uq plays the role of the number of particles in the Wigner-Seitz cell. 
Under the above conditions, seen as a function of the ^-variable, the bulk zero-field orbital sus- 
ceptibility at fixed /3 > and density of particles po{R), i? > is defined by: 

XniP, po{R),b - 0) := Xr{P, po{R),b = 0), 6 = 0). (1.16) 

Before giving our main result which is concerned with the quantity defined in (|1.16p . let us 
introduce some more notation. On C^(R"^) define the 'single atom' Schrodinger operator: 

Hp:=^i-iVf+u, (1.17) 

where u is the function which appears in (jl.7p and obeys assumption (i/r)- Then Hp is essentially 
self-adjoint and its self-adjoint extension (denoted again by Hp) with domain ■H^(M^), is bounded 
from below. Furthermore, acssiHp) = [0, oo) is absolutely continuous, and Hp has finitely many 
eigenvalues in (— oo,0) if any, see e.g. [HI Thm. XIII. 15]. Throughout this study, we suppose: 

(jz/ni) Hp has at least one eigenvalue in (— oo,0), 

together with the non-degeneracy assumption: 

(•^nd) All the eigenvalues of Hp in (— oo,0) are simple (i.e. non- degenerate). 
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Then wc denote by {Xi}J^i, r G N* the set of eigenvalues of Hp in (—00, 0) counting in increasing 
order, and by {^i}J^i the set of corresponding normahzed eigenfunctions. As weU, we denote by 
Hi the orthogonal projection onto the eigenvector and we define 11^ :— 1 — 11;. 

In the presence of the uniform magnetic field (as in (|1.9p ). define on C5"(R'^) the 'single atom' 
magnetic Schrodinger operator: 

Hp{b) ■.= ^(-i\/ -baf +u, beR. (1.18) 

By [171 Thm. B.13.4], V6 G M p.lSp is essentially self-adjoint and its self-adjoint extension, 
denoted again by Hp{b), is bounded from below. The nature of the spectrum of Hp(b) is not 
known in general, except for b small enough. Indeed from [SJ Thm. 6.1], the eigenvalues of Hp in 
(—00, 0) are stable under the perturbation by the magnetic field provided that it is weak. From 
the asymptotic perturbation theory in [35l Sec. VIII] and due to the assumption {s^^d), then there 
exists b > s.t. V|&| < b and any / G {1, . . . , r}, Hp{b) has exactly one and only one eigenvalue 
\i{b) near A; which reduces to A; in the limit & — 0. In particular, each eigenvalue Xi{- ) can be 
written in terms of an asymptotic power series in &, see e.g. [71 Thm. 1.2]. Hereafter we denote 
by {Ai(6)}[^j^, T G N* and |6| < b, the set of these eigenvalues for Hp{b) counting in increasing 
order. 



We now formulate our main result. By [571 Thm. B.7.2], for any ^ G C \ [\aia{Hp), c«), the 
resolvent operator {Hp — ^)~^ is an integral operator with integral kernel {Hp — ^)~^(- , • ) jointly 
continuous on M*^ \ L), £) := {(x, y) G M*" : x = y}. Introduce on L'^{M?) the operators Tpj{£), 
j = 1,2 generated via their kernel respectively defined as: 

V(x,y)GM^\A rp,i(x,y;0 :=a(x-y).(*Vx)(ffp-0~'(x,y), (1.19) 
rp.2(x,y;0 :-ia2(x-y)(ifp-0-'(x,y). (1.20) 

Our main result provides an asymptotic expansion in the tight-binding situation (i.e. when 
the distance R between two consecutive ions is large) of the bulk zero-field orbital susceptibility 
defined in (jl.l6p when the number of particles in D,p is fixed and in the zero-temperature limit: 

Theorem 1.1 Suppose {s^t), {■^m) cmd (^nd)- Assume that the number of particles tiq G N* in 
the Wigner-Seitz cell is fixed and satisfies uq < t, while the density is given by (|1.15p . Then: 
(i). For any < a < 1 there exists a R-independent constant c > s.t. 

Xr{po{R)):= limA'p(/3,po(i?),6 = 0) = -^jrp(no) + 0(c-'=^°), (1.21) 

Ptco |SZp| 

with: 

2 . 



^p(no) 



^) lTri.(K3)|^ d^aHp-0-'[TpA{0Tp,i{0-Tp,2{0]'j, (1-22) 



where T„g is any positively oriented simple closed contour enclosing the uq smallest eigenvalues of 
Hp while letting the rest of the spectrum outside. 

(ii) . The R-independent quantity in the leading term of the expansion (|1.2ip can be identified with: 

\ 2 no ,2 , 

q\ d^A; 



.rp{no) = -m E^(^ = 0)- (1-23) 



1=1 



(iii). The leading term of the expansion (|1.2ip can be rewritten as a sum of two contributions: 

^ ^p(no) = ^^La(no) + ^^yy(no), (1.24) 
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with: 



^Laino):^~(^] -^^($;,(X2+X2)<I,,), (1.25) 



j^^^vvM - (y ^^(Xa^.ln-lHp-AOn-l-La*.). (1.26) 

_ffere := X • 6^. k e {1,2,3} stands for the position operator projected in the k-th direction, 
and Lfc := L • ek the k-th component of the orbital angular momentum operator L := X x (— iV). 

Remark 1.2 The contribution in (|1.25p is the so-called Larmor diamagnetic susceptibility. It is 
generated by the quadratic term of the Zeeman Hamiltonian. Still assuming one single ion in the 
Wigner-Seitz cell, it reduces in the classical limit to the Langevin formula, see (jl.ip ; 

\ 2 -. n o 

where ri denotes the distance from the origin to the l-th particle in the plane orthogonal to 63 (as a 
rule the center of mass of the ion nucleus is located at the origin and the fixed axis passing through 
it is taken parallel to the magnetic field). As for the contribution in (jl.26p . it is the so-called 
Van Vleck orbital susceptibility. It involves the third component of the orbital angular momentum 
operator which comes from the Zeeman interaction energy. We mention that if uq ~ 1, then the 
orbital Van- Vleck susceptibility is purely paramagnetic. Otherwise, (jl.26p can be written as a sum 
of a positive and negative contribution. 



1.3 Discussion on the assumptions. 

Let us first discuss the assumption (is^J ) on the single-site potential u. The physical modeling 
requires that the 'single atom' operator Hp in p.l7p possesses finitely many eigenvalues (at least 
one) below the essential spectrum; therefore u has to be chosen accordingly. Choosing u compactly 
supported plays an important role in our analysis. This latter is based on an approximation of the 
resolvent operator {Hr — C)^^ ii^ the tight-binding situation via a geometric perturbation theory, 
see Sec. [Hand (|2.13p . The fact that u has a compact support leads to exponentially decreasing 
estimates in R°' , < a < 1 on the 'error term' arising from the approximation of {Hp — ^)~^ with 
{Hp — ^)^^ in the bulk of the Wigner-Seitz cell for large values of R; see discussion below (|2.13p . 
This gives rise to the exponentially decreasing behavior in i?" of the remainder in the asymptotic 
expansion ()1.21|) . see Proposition 13.71 However, we believe that the optimal a's is a = 1, i.e. the 
remainder should behave like 0{e~'^^), c > 0. We think that it could be obtained from our analysis 
by using a more refined geometric perturbation theory to approximate the resolvent {Hp — ■ 

Furthermore, we point out that the leading term in the asymptotic expansion (|1.2ip is un- 
changed if the single-site potential is not chosen compactly supported. Consider the assumption: 

{£!/,.) u e C\R^;M.) with u = 0{\x\-^^+^'^) for |x| sufficiently large. 

Let us note that from [441 Thm. XIII. 6], acss{Hp) = [0, 00) and Hp has a finite number of bound 
states in (— oo,0). Replacing (jz/r) with (s^r') in the assumptions of Theorem ll.il then under the 
same conditions, one may expect the behavior of the remainder to be only polynomially decreasing 
with R owing to the 'tail' of the potential. Finally, let us mention that u is chosen continuously 
differentiable to ensure some regularities for the eigenvectors of Hp , see (|4.15p and Lemma 13.141 

Let us discuss the non-degeneracy assumption (i/nd)- Our analysis is based on the insulating 
situation which occurs when the Fermi energy lies in the middle of a spectral gap of Hp, see [111 
Thm 1.1]. When the number of particles ng in the Wigner-Seitz cell is any integer lesser than 
the number of negative eigenvalues of Hp while the density is given by (|1.15p . this together with 
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(^nd) automatically lead to the insulating situation for R sufficiently large, see Proposition 13. II 

Nonetheless, we stress the point that when getting rid of the assumption (^nd), then the 
insulating condition can still occur for R sufficiently large provided that one sets some restrictions 
on the rio's in (jl.lSp . It has to obey (see proof of Proposition 13.11 and Remark l3.5p : 

no < T and 3>c G {1, . . . , i^} s.t. no = dim (oi . (1-28) 

1=1 

Here r is the number of all the eigenvalues of Hp in (— oo,0) counting multiplicities, < r is the 
number of distinct eigenvalues and (oi, I E {l,...,z^} stands for the eigenspace associated with 
the (possibly degenerate) eigenvalue A; oi Hp, Xi = {A[™^}^™f'. Supposing only (^■)-(^m), and 
assuming that the number of particles uq in the Wigner-Seitz cell is fixed and obeys (jl.28p . while 
the density is given by (|1.15|) . then the formulae in Theorem 11.11 (ii)-(iii) have to me modified 
accordingly (the statement in (i) is unchanged). Thus (|1.23p becomes: 

/ ^ 2 K dirndl ,2 \ (™) 
^ 1=1 m=l 

and (|1.25p - (|1.26p respectively become (with intuitive notations): 

/ X 2 , J< dimSi 

\^\^vvM--- (?) E (i3<i>|'"\{n|'")'^(Fp-A|'"))n|™)'^}-iL3<i>|")). 



In a way, the condition in (|1.28p can be linked with the one concerning the filling of the electron 
shells mentioned in Sec. II. II when dealing with the susceptibility of an ion, see below (|1.6p . When 
all its electron shells are fulfilled, it is found that the zero-field orbital susceptibility reduces to 
the diamagnetic Larmor contribution. But this result leans on the Hund's rules which state that 
such an ion has necessarily its total electronic orbital angular momentum null in its ground state. 
Putting aside these considerations of a purely atomic nature, then the Van-Vleck contribution has 
to be taken into account. 

1.4 An open problem. 

From the foregoing, a natural question arises: does the insulator condition occur when no is 
any integer less than the number of negative eigenvalues of Hp, counting multiplicities if some 
of the eigenvalues are degenerate, while the density is given by (|1.15p ? Such a problem comes 
up when the single-site potential is chosen spherically symmetric. To tackle it we need to know 
precisely the behavior of the negative spectral bands of Hp near the degenerate eigenvalues of 
Hp in (— oo, 0) for large values of R. For instance suppose that one of the negative eigenvalue of 
Hp, say Ac, is two-fold degenerate. For R sufficiently large, it is well-known that the spectrum 
of Hp in a neighborhood of Ac consists of the union of two Bloch bands, see Lemma and also 
[211 Thm. 2.1]. But we need to know much more; in particular we need to control how the Bloch 
bands behave the one relative to the other. Especially, do they always overlap for R sufficiently 
large or do they overlap only in the limit R f oo? How fast each of the Bloch bands reduce to the 
Ac's? This remains a challenging spectral problem. 

1.5 The content of the paper. 

Our current paper is organized as follows. In Sec. [2l we use a geometric perturbation theory 
to approximate the resolvent operator {Hp — the tight-binding situation. The key idea 
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consists in isolating in the Wigner-Seitz cell the region close to the boundary from the bulk where 
only the 'single atom' operator in (|1.17p will act. The use of well-chosen cut-off functions allows 
us to keep a good control on the 'error term' arising from the approximation of {H^ — ^)~^ with 
{Hp — C)^^ in the bulk of This will play an important role to obtain the behavior with R 
of the remainder in the expansion (|1.2ip . Sec. [3] is devoted to the proof of Theorem 11.11 In 
Sec. 13.11 we show that under our assumptions only the insulating situation can occur in the 
tight-binding situation. It is based on the Bloch-Floquet theory and the behavior of the negative 
spectral bands of Hp near the negative eigenvalues of Hp for large values of R. Subsequently, 
we prove the asymptotic expansion in (|1.2ip . The method requires the insulating situation so 
as to perform the zero-temperature limit in the involved quantity. The main ingredients are the 
approximation from Sec. [2] and the exponential localization of the eigenvectors associated with 
the negative eigenvalues of Hp. In Sec. 13.21 we prove the identity (|1.23p . It is based on the 
stability of the negative eigenvalues of the operator Hp when considering weak magnetic fields. 
The crucial ingredient is the gauge invariant magnetic perturbation theory, applied to the kernel of 
the resolvent operator of Hp{b), allowing to keep a good control on the linear growth arising from 
the vector potential. In Sec. I3.3l we prove the formula (jl.24p . It is based on the so-called Feshbach 
formula and requires the use of the asymptotic perturbation theory for the singular perturbation 
Hp{h) — Hp to control the behavior of the involved quantities for small values of h. In Sec. HI we 
have gathered together all the proves of the technical intermediary results needed in Sec. [31 

2 An approximation of the resolvent via a geometric per- 
turbation theory. 

The method we use below is borrowed from [ini [TT] . 

For any < a < 1, < K < 3 and i? > define Qr{k) {x e IT^ : dist(x, Sf^i?) < Ki?"}. 
Then for R sufficiently large, models a 'thin' compact subset of Vtp near the boundary with 

Lebesgue-measure of order 0(i?^+°). 

Let < a < 1 be fixed. Below by R sufficiently large we mean R> Rq with Rq = Ro{a) > 1 s.t. 



(2.1) 



Let us introduce some wcU-chosen family of smooth cutoff functions. 
Let gp and gp, R > Rq satisfying: 



Snppigp) c{nR\eR{l)), 0<gR<l; 
Supp(gi?J C (f7fl\ 9^(1/2)), = 1 when xe (r!fl\ 9^(1)), < .g^?, < 1. 



Moreover there exists a constant C > s.t. 
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Vi? > i?o, max- 



V|5|<2, |5|=^S,. 



i=i 



With these properties, one straightforwardly gets: 



dist(Supp(i?^5fl), Snppigp)) > CR"^, VI < |s| < 2, 



gngR = 9ri 



(2.2) 
(2.3) 



for another i?-independent constant C > 0. Also let gp, R> Rq satisfying: 



Supp(lfl) C {np \ 9«(3/2))C, Ip = 1 when x e (17^ \ 9«(1)) 



c 



< 5fl < 1; 
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for another _R-indepcndcnt constant C > 0. Note that with these properties: 

^9R{l-gB) = a-9R), (2.4) 
dist(Supp(i:>"|^), Supp(l - qr)) > CR°', VI < \s\ < 2. (2.5) 

Let us now define a series of operators. At first, introduce Vi? > Rq and V$ £ q{Hr) n g{Hp) 
(here g{- ) denotes the resolvent set) on L^(M'^): 

7^i^(0 9RiHp ^ ^)-'9R + hiHR ~ - gn). (2.6) 

In virtue of the support of cutoff functions, one has: 

{Hr - OijR = {Hp - 09R- 

Using that Ran(7?.i?(^)) C Dom(iJj^) by standard arguments, then from (|2.2p along with (|2.4p : 

{HR-onpiO^i + WRiO, 

where, Vi? > Rq and G giHp) n ^(iJp): 
Wii(e) := {-^(A5i?) - (V5h) • V}(i7p - 

+ {-i(A|p) - {V§r) ■ V}{Hr - e)"'(l - 9r)- 
Since Wi?,(^) is bounded, see e.g. [12l Lem. 5.1], this means that in the bounded operators sense: 

{Hr - = 7^i^,(C) - {Hr - 0-'Wr{0- (2.7) 
Afterwards, 72,ij(^) in p.6p can be rewritten via the use of the second resolvent equation: 

nR{O^^R{0-^R{0, (2.8) 

where, Vi? > Rq and V^ G g{HR) n g{Hp): 

^r{0 ■■= <3r{Hp - 0~'gR+h{Hp - e)-'(l - 9r), (2.9) 

'^r{0 ■■= h.{HR - 0''Vr{Hp - - 9r), (2.10) 

with: 

Vr := - ^^)- (2.11) 

uGZ3\{0} 

Finally ^r{£,) in (|2.9p can be rewritten: 

^r{0^{Hp-0'^+'^r{0, (2.12) 
where, Vi? > i?o and V^ € £i(iip) n ^(iip): 

2Bp(0 := {Hp - e)-M-^(Agp) - {Vcjr) ■ V}{Hp - 0'^gR+ 

+ {Hp - e)-M-^(A|p) - (V|p) • V}{Hp - e)-'(l - 9r)- 
Gathering (|2.7p . (|2.8p and (|2.12l) together, we get in the bounded operators sense on L^(R^): 

{Hr - - {Hp - 0^' + 2Bp(0 - ^r{0 - {Hr - 0''Wr{0- (2.13) 
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Let us now justify why the first term in the l.h.s. of (|2.13p turns out to be a suitable approxi- 
mation when R is sufficiently large. To do that define i?i > 1 so that: 

Snppiu) C (ilnA^mi^))- (2.14) 

A such i?i's exists since the support of u is compact, see assumption {^/i). Now look at the r.h.s. 
of (|2.13p . At first, Vi? > Rq the last term {Hn — 0~^^RiO h^-^ its operator norm exponentially 
small in i?", see p.23p below. The same holds true for the operator norm of Wr{^), see ()2.23|) too. 
However, this is not the case for the operator norm of #^(0 even for R> Rq and large enough, 
see (|2.22p . This comes from the fact that the support of (or (1 — gji)) and the support of Vr 
in (|2.1ip are not disjoint. But VR > max{i?o,-Ri} the Hilbert-Schmidt norm of XiiR^niO ^^.d 
^niOxnR, where Xfin denotes the indicator function of fljf,, are exponentially small in R", see 
(|2.24[) . This last feature, which results from the fact that u is compactly supported, will turn out 
to be decisive to get the exponential decay of the remainder in the asymptotic expansion p.2ip . 
We end this paragraph by giving a series of estimates which we will use throughout: 

Lemma 2.1 Let S = P or R. For every rj > (and Vi? > when 'E. = R) there exists a constant 
•& = i?(?7) > and a polynomial p{- ) s.t. G C satisfying dist(^, CT(i?H)) > i]: 

V(x,y)eR«\A |(gH-e)-^(x,y)|<p(iei) V''^T' ' ^^---TirWv (^.IS) 

|x-y| i + Kl 

|Vx(i/H - 0"'(x,y)| < pmf-^^- (2.16) 

Proof. (Pl^ and ^^^Y^ follow from gll Thm. B.7.2] and [H Lem. 2.4] respectively. □ 

Lemma 2.2 Let < a < 1 he fixed and i?o — Ro{a) > 1 as in p.ip . Then for every ij > there 
exists a constant = "d^rj) > and a polynomial p{- ) s.t. 

(i) . Vi? > Ro and G C satisfying dist(^, cr(iJfl) n a{Hp)) > rj: 

V(x,y)eR6\A max{|(7^^(C))(x,y)|,|(^^(0)(x,y)|}<p(|e|)^^^i^, (2.17) 

|x-y| 

max{|Vx(7^«(0)(x,y)|,|V.(,^«(0)(x,y)|}<p(|^|)^^^^. (2.18) 

(ii) . Vi? > i?o andyS, e C satisfying dist(^, cr(iiii) n cr{Hp)) > 77; 

V(x,y)eM«, |(^«(0)(x,y)| <p(|CI)e-''«l''-''l, (2.19) 
max{|(Wfi(0)(x,y)|,|(2nfl(e))(x,y)|,|Vx(2BH(C))(x,y)|}<p(|C|)e-''«^"e-^^l''-'^l. (2.20) 

(iii) . Vi? > max{i?o, i?i} (see ((^1^ ) and V$ € C satisfying dist(C, (j{Hr) n a{Hp)) > rj: 

V(x,y) G R^ max{|(xo„ii?j(x)(ii;j, - 0"'(x, y)i>i?(y)|, 

|FH(x)(iip - e)-^x,y)(xo Jl - .gi^))(y)|} < P(l^l)e-''^^"c-''«l--''l. 

Proof. Let us prove (i). (|2.17p directly follows from (|2.15p . On the other hand: 

Vx(7^i^,(0)(x,y) = {(V,te)(x)(iip -0"'(x,y) +gi?(x)Vx(iip -0"'(x,y)},te(y)+ 

+ {(V|p)(x)(iip - 0"'(x,y) + |i?(x)Vx(ifp - 0"'(x,y)}(l - gniy)), 

and by replacing (Hp - ^)-i with {Hp - ^)-^ above, we get Vx(.^p(0)(- , ' ). Then ^^JE^ follows 
from the properties (|2.3p - p.5p together with the estimates (|2.15p - (|2.16p . Let us prove (ii). The 
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first estimate results from tlie assumption (^) and (|2.15p . see [TH Lem. A. 2]. Due to (|2.3p - (|2.5p 
again together with (j2.15p . then under the conditions of the lemma (below S = P or i?): 

max{ \{Agn) (x) {H^ - C) " ^ (x, y)^^ (y ) | , | ( V^h) (x) {H^ - 0^' (x, y)9R (y) | , 

|(A|H)(x)(i/H - C)-'(x,y)(l - gR){y)\, |(v|«)(x)Vx(ifH - C)-'(x,y)(l - .gi^)(y)|} 

<p(|^|)e-''«^°e-''«l''-yl, (2.21) 

for another i?-independent t? > and polynomial p(-). This leads to (|2.20p for the kernels of 
Wfl(C) and 2IJi?(0- As for Vx(22Jfl(0)(- , • )> it is sufficient to use ([2211), (Pl^ along with [H 
Eq. (A.12)]. Finally (in) follows from ((2TT4)) ensuring that dist(Supp(u), Supp(xn^|_R)) > ^R^.n 

Remark 2.3 (i). From Lemma \2.1\ and Lemma \2.2\ (i)-(ii), together with the Shur-Holmgren 
criterion, one has \IR > Ro : 

max{||(i7fl^ - O-'ll, ||7^fl,(OIUI^fl(^)IU|V7^fl(OIU|V^fl(0||, ll^flXOII} < P(iei), (2.22) 
max{|lWfl(0IU|2ITfl,(e)|| || V2IT;^,(e)l|} < p(kl)e-''«''° , (2-23) 
for another R-independent constant ■& > and polynomial p( - ). 

(ii). Let (32(i^(R'^)), II ■ II32) ^^'^ Banach space of Hilbert- Schmidt operators. By using the 
* -ideal property of 32{L'^ {S?)) , then from Lemma \2.S\ (iii) one has Vi? > max{i?,Oi ^i}-' 

uis^^{\\xnMm-j.A\^RiOxnJ,A < MI^De"''^''", (2-24) 
for another R-independent constant 1!} > and polynomial p{- ). 



3 Proof of Theorem 11.11 

This section is organized as follows. The first part is devoted to the proof of the asymptotic 
expansion ()1.2ip in the tight-binding situation and in the zero-temperature limit. The second and 
third part are respectively concerned with the proof of the identities (jl.23p and (jl.24p . For reader's 
convenience, the proof of technical intermediary results are collected in Appendix, see Sec. |4l 



3.1 Proof of (i). 

3.1.1 The location of the Fermi energy. 

Here we are interested in the location of the Fermi energy in the tight-binding situation when 
the number of particles no G N* in the Wigner-Setz is fixed, while the density is given by (|1.15p . 
Recall that under our conditions, Vi? > the Fermi energy: 

£R,FipoiR)) ■■= lim Mfl'(/3,Po(i?),& = 0), 
always exists, see [TTl Thm 1.1]. 



Before giving the main result of this paragraph, let us introduce some notations within the 
framework of the Bloch-Floquet theory. For details, we refer to [6l Sec. 3.5] and [52]. The results 
we give below hold true for any R > 0. Denote by fi^ the unit cell of the dual lattice {2-k/R)'I? 
(the so-called first Brillouin zone) of the Bravais-lattice Rl? . With 5(]R'^) denoting the Schwartz 
space of rapidly decreasing functions on R^, consider the Bloeh-Floquet(-Zak) transformation: 

U -.SiW^)^ L'^{VL}^,L'^{nR))'^ [ dkL2(r2fl) 

(Z^0)(x;k) = ^l= ^ e-'''-(^+-)0(x + t^), k e r!Jj, X e f^fl^, </) G 5(M3), 
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which can be continued in a unitary operator on L^(M^). The unitary transformation of Hn is 
decomposable into a direct integral lAHnU* ~ J^, dk/ifl(k), where: 

lives in L^(M^/i?Z'^). By standard arguments, hji is essentially self-adjoint on C°°{M.'^ / RZ^); the 
domain of its closure is the Sobolev space H'^ {M.^ / RZ^) . For each k g /ifl(k) has purely 
discrete spectrum with an accumulation point at infinity. Then we denote by {£'/?.((k)};>i the 
set of eigenvalues counting multiplicities and in increasing order. Due to this choice of labeling, 
the Ejij^s are periodic and Lipschitz continuous on f2Jj. Indeed they are not differentiable on a 
zero Lebesgue- measure subset of fijj corresponding to crossing-points. If Z > 1, the l-th Bloch 
band function is defined by Erj := [minkenjj (k), maxkeoj, ^'^.((k)]. The spectrum of Hr is 
absolutely continuous and given (as a set of points) by cr{Hji) = Ui^i'^^?,'- Note that the sets 
Sfi^i can overlap each other in many ways, and some of them can even coincide. The energy bands 
are disjoint unions of Sr/s. Moreover, if max Sr^i < minf/j^j+i for some I > 1 then we have a 
spectral gap. Since the Bethe-Sommerfeld conjecture holds true under our conditions, see e.g. pSl 
Coro. 2.3], then the number of spectral gaps is finite, if not zero. 

Our main result below states that, under our conditions and in the tight-binding situation, the 
Fermi energy always lies in the middle of a spectral gap of Hr (i.e. only the insulating situation 
can occur), and moreover, it provides an asymptotic expansion of the Fermi energy; 

Proposition 3.1 Suppose (^), (Mn) and (^id)- Assume that the number of particles hq G N* 
in the Wigner-Seitz cell is fixed and satisfies uq < t, while the density is given by (jl.lSp . Then: 

(i) . For any {3 > 0, let fj,^^\/3, p(j{R),b = 0) G M 6e the unique solution of the equation pR{l3,p,b = 
0) = Pq (R) . Then for R sufficiently large, the Fermi energy satisfies: 

SRApoiR)) lim pf{P,poiR),b = 0) = +"^i^g^,»o+i < , 

pjoo Z 

(ii) . Let us define: 

P \ ■ / (^«o +^«o+i)/2 when no <T, 

£p,fM I ;^^/2 when no = r. ^^■^) 

Then under the additional assumption that uq < t, one has: 

£rApo{R)) = £p.F{no) + ©(c-Vl^-o+ilfl). 

Remark 3.2 (i). We stress the point that the non- degeneracy assumption (.s^i^d) along with the 
fact that po{R) is given as in (|1.15p . together imply the insulating situation for R sufficiently large. 
(ii) . We will see in Proposition lS. that the Fermi energy actually plays any role for the statements 
of Theorem \l.l \ since it can be remove of the main quantities without changing their values. 

The rest of this paragraph is devoted to the proof of Proposition 13.11 
Let us start by writing down an expression for the bulk density of particles. Under the grand- 
canonical conditions, let /? > and /i £ M. For any i? > 0, let '^jf'^ be the counter-clockwise 
oriented simple contour around the interval [inf (t{Hr),oo) defined by: 

:= m € ['5fl,oo), 3^ = ±^}U{5Re = -5^, 5^ e [-^, ^]}, 5r := \nia{HR) - 1. (3.3) 

Let us note that for any R > Q, the closed subset surrounding by is a strict subset of the 

holomorphic domain S := {C G C : 3C e (— 7r//3, 7r//3)} of the Fermi-Dirac distribution function 
^FD{P,pi]£.) := e'^('^-«)(l -He^(^-«))-^. From (11131) and seen as a function of the ^-variable, the 
bulk zero-field density of particles reads V/3 > 0, G M and Vi? > as, see e.g. [12l Eq. (6.3)]: 

Pfl(/3, 6 = 0) := ^Tri,(^3) jxn,. j^^^^ d^ ^fd{P, P; 0{Hr - 0"') | • (3.4) 
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Wc mention that another way to express the bulk zero-field density consists in bringing into play 
the integrated density of states of the operator Hr. Under the conditions of (|3.4p . 

/i, 6 = 0) = - dt ^{t)NR{t), 

where Nr{- ) denotes here the integrated density of states of Hr = HrQ) = 0) defined in (jl.lOp . 
We recall that when the magnetic field vanishes, Nr is a positive, continuous and non-decreasing 
function, and it is piecewise constant when the energy parameter belongs to a spectral gap. 
In order to write down an expression for the bulk zero- field density in the zero-temperature limit, 
we need to rewrite p.4p by the use of the Bloch-Floquet decomposition, see e.g. [HI Sec. 2]. In 
view of our notations introduced above, we collect in the following lemma all the needed results: 

Lemma 3.3 (i). Let (3 > Q and ^ G M. Then for any R > 0; 

1 °° /■ 

PR{l3,fi,b = 0)^ ||o*| E / , dkfFDW,^^;ERAk)). 

I -R N R\ j — 1 

(ii). For any R > 0, let ii > inf (j{Hr) he fixed. We have the identity: 

lim pR{P,fi,b:^ 0) = V / dkxiinf aiHn).f,]{ERj{k)) = Nr{p), (3.5) 

where X[in! (j{Hii)4i]{') denotes the indicator function of the compact interval [ini a{HR), p] and Nr 
the integrated density of states of the operator Hr . 

Now let us get back to the location of the Fermi energy in the tight-binding situation when the 
density is given by (|1.15p . We need first to know how the negative spectral bands of the operator 
Hr are localized at negative eigenvalues of the operator Hp for large values of the i?-parameter. 
For the sake of completeness, in the below lemma we do not consider assumption (^nd) and then 
we allow the negative eigenvalues of Hp in (— oo, 0) to have some degeneracies: 

Lemma 3.4 Let r G N* &e the number of the eigenvalues of Hp in (— oo, 0) counting multiplicities. 
Denote by {A(}^^j^ with v < t the set of distinct eigenvalues counting in increasing order. Then 
for R sufficiently large, there exist real numbers — Ai ^ \ < cr^i, dp^i < 0, Z = 1, . . . , z/ and Cr^^+i 
satisfying cr^i < dps < ■ ■ ■ < cr^^ < dp^^ < Cp^^+i s.t. 

(i) cr(i?i?)r(_oo,0) ^\fi=i[cR,udRj], 
(u) [cR^i,dRj] n [cfl,„i, d_R,m] ^% for I m, 

(iii) A; e icR,i,dRj), 

(iv) dp^i, + Cp^y+i < 0, 

together with the properties that cpj, dp.i — > A;, I = 1, . . . ,v and Cp^^+i — >■ when R^ oo. 
Furthermore we have the following relation with the Block bands of Hp : 

[cRj,dpj]n a{Hp) = (J Ep^rn, l^l,...,iy, 

7H — 1 

where Si stands for the eigenspace associated with the possibly degenerate eigenvalue A; . 
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Proof. These statements directly follow from Thm. 2.1] taking into account our choice of 
labeling for the Eu/s (i.e. increasing order). Note that Theorem 2.1 in ^T] is established under 
the assumptions that Vr is smooth and sufficiently fast decaying at infinity. But the statements 
still hold true under our conditions on Vr, see [171 Thm. 2]. □ 



Now we are ready to prove Proposition 13. II 
Proof of Proposition [37ll Let us first prove (i). Consider the equation: 

Trrrryi / dkx[inf<T(if«),iJ](-Bi?,i(k)) = no. 

Due to the non-degeneracy assumption {^/nd), Lemma [3^ ensures that the Bloch bands £r^i, I = 
1, . . . , T are simple, isolated from each other and from the rest of the spectrum for large values of R. 
Hence if rig e N* satisfies tiq < t, then E must belong to [maxkgnj, -E'fl.noC*), minkgnj, -E'fl_„j,+i(k)]. 
This comes from the fact that the Lebesgue-measure of the set {k e iljj : -E;(k) < E} equals |fi^| 
if and only if i? > maxkGiijj -^((k). Getting back to p.5|) . this means that for large values of R: 

po{R) ^ Nr{E) V-B G [max£_R,„o,min£_R,„o+i]. (3.6) 

Now p.ip follows from [11] Thm. 1.1]. Let us prove (ii). Supposing that no < r, use that: 

maxfij_„o — A„o min£ij^„o+i — A„j,+i 



■'R,F 



2 2 

together with the following estimate which holds uniformly in k S see e.g. [Ul Thm. 2] 



\ERj{-k)\ - = 0(i?-ie-Vl^n^), / = 1, . . . ,r. 

□ 

Remark 3.5 Let us give a sufficient condition ensuring the insulator condition when considering 
some degeneracies for the eigenvalues of Hp in (— oo,0). From Lemma \3.4\ then for R sufficiently 
large (|3.6p holds Vi? belonging to a spectral gap of Hr provided that: 



3>f G {1, . . . , I'} s.t. Tip = dim (g; . 



1=1 

3.1.2 Isolating the main contribution at zero-temperature. 

We start by writing down an expression for the bulk zero-field orbital susceptibility. Under 
the grand-canonical conditions, let /? > and p e M. For any i? > let be the positively 

oriented simple contour around the interval [inf a {H r) , oo) defined in (|3.3[) . Then Vi? > the 
closed subset surrounding by '^p^^' is a strict subset of the holomorphic domain S := {C e C : 
5CG (-7r//3,7r//3)} of the map ^ f(/3,^;^) := ln(l + e'^^^"?)). Let us note that f(/3, /x; • ) admits 
an exponential decay on i.e. there exists a /3-independcnt constant c > s.t. 

V^e^^f\ |f(/3,M;C)l<ce'''^e-^'^«. (3.7) 

From (|1.14p and seen as a function of the ^-variable, the bulk zero-field orbital susceptibility reads 
V/? > 0, V// e M and Vi? > as, see e.g. [Ill Eq. (1.21)]: 

XR{j3,^l,h^Q)■.^ f^) / def(/3,/i;0x 

\cj l3\ilR\2T: Jr^w 

X Tri.(K3){xn„(iii? - S.r^[TR,i{i)TRAi) - Th,2(C)]xoh}, (3.8) 
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where Tfl j(^), j — 1,2 arc bounded operators generated via their kernel respectively defined as: 

V(x,y) eM*^\Z?, rfl,i(x,y;0 :=a(x-y).(iVx)(Hfl-e)"'(x,y), (3.9) 

rfl,2(x,y;0 :=^a2(x-y)(iffl,-0"'(x,y). (3.10) 

Recall that we have introduced the operators Tp,j{£^), j = 1,2 via their kernel defined similarly to 
(I311)-(ISIini)butwith(ffp-C)-Mnsteadof(i/i?-0-\sec([rini)-([r2n^ Since |a(x-y)| < |x-y|, 
then under the conditions of Lemma [2TT] (below S := i? or P): 

yi^,y)eR'\D, |r.^.(x,y;OI <p(|e|)^^^, := j = 1,2, (3.11) 

for another constant d > and polynomial p{- ). Due to the estimates p.lip and p.lSp . the op- 
erators {He - 0~^Ts,iiOTs,i{0 and [Hs - ^)"^7h,2(C), S = i? or P both are locally trace class 
on L^(R'^). Furthermore, both are integral operators with a jointly continuous integral kernel on 
MP, whose diagonal part is bounded above by some polynomial in |^| uniformly in the spacial vari- 
able, see e.g. [T^ Lem. A.l]. This along with p.7p ensure that the quantity in p.8p is well-defined. 

Now let us turn to the actual proof of p.2ip . We point out that the main difficulty consists 
in isolating the main contribution from p.8p in the tight-binding situation, while keeping a good 
control on the behavior of the 'remainder' term, even in the zero-temperature limit. 

The starting point is the approximation of the resolvent operator {Hji — S,)^^ derived in Sec. 
[2l By replacing in (|3.8|) each resolvent {Hji — ^)~^ (look at the definitions (|3.9p - p.l0p ) with the 
r.h.s. of p.lSp . and taking into account the features of the three last operators of the r.h.s. of 
(|2.13p we discussed in Sec. [21 we naturally expect the main (still P-dependent) contribution from 
p.Sp to be obtained by replacing each operator {Hji — ^)~^ with {Hp — ^)~^. In this way define, 

> 0, V/x e R and VP > the following quantities, 

Xn{P,^^,b^O)■.^ (^) i-i^- /^^ def(/3,Ai;e)x 

X Tri.(K3){xojPp - O^'ITpAOTpAO ~ Tp^2mxn^}, (3-12) 

ApiP, fi) := XpiP, M, & = 0) - Xr{P, fi,b^O), (3.13) 

where ^j^"* in (|3.12p denotes the contour defined as in (|3.3p but with Sp := mf(a{Hp)) — 1 = Ai — 1 
instead of Sp. 

Now consider the canonical conditions. Suppose that the number of particles in the Wigner- 
Seitz cell rig € N* is fixed and obeys no < r, while the density is given by (|1.15p . Let /i^^ (/3, po{R), b = 
0) G K be the unique solution of the equation pR{f3, e^^, 6 = 0)= Pq[R). Then from (|1.16p together 
with ((3T2l) - ((3T3| . one has > and VP > 0: 

Xr{P, poiR), 6 = 0)= Xp{l3, pf {13, po{R), 6 = 0), 6 = 0) + Ap{p, pf (/3, po{R), 6 = 0)). (3.14) 

The next step of the proof consists in performing the zero-temperature limit in ()3.14p in the tight- 
binding situation. Here the crucial point is the insulator situation: for P sufficiently large, the 
Fermi energy lies outside the spectrum of Hp and is located in a neighborhood of the middle of 
the interval (A„o,A„q+i) if no < r, (At-,0) otherwise; see Proposition 13. II along with Lemma [3.41 
Remind that the insulator condition results from the non-degeneracy assumption {d^nd) together 
with the condition (jl.lSp . To perform the zero-temperature limit in p.l4p . we need the two 
following results whose proves can be found in Appendix, see Sec. 14.11 We recall that {Xi}J^i, 
r G N* stands for the set of eigenvalues of Hp in (— oo, 0) counting in increasing order. 
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Proposition 3.6 Let /<;, <; G {1,...,t} be an open interval s.t. C (A;;,A;;+i) and {\ + 

\j^i)/2 € when <; < t; otherwise Ir C (Ar,0) and Ar/2 G /r- T/ien Vi? > and for any 
compact subset K C I^: 

lim XR{P,t,,b^ 0) ^ -^XRi^i,b^ 0), 

/3too \Ur\ 

uniformly in fi € K , with: 

<?j^(/i,6 = 0) (^^^ de(/.-e)x 

where stands for any positively oriented simple closed contour surrounding the smallest eigen- 
values of Hp in (—00,0) while letting outside the rest of the spectrum of Hp. 

Proposition 3.7 Let /<;, G {1, . . . , r} &e an open interval as above. Then for R sufficiently large, 
lim^-f-oo Aij(/3, /^) exists uniformly on compact subsets K (Z L^. Denote it by Aij(/i). Furthermore 
VO < a < 1 there exist two constants c,C > s.t. V/.t G /<, and for R sufficiently large: 

|Ap(A.)| <C(l + H)e-=«°. (3.15) 

Let us emphasize that the exponentially decaying estimate appearing in p.lSp arises from the 
fact that u is compactly supported, see assumption (^r)- The proof of Proposition 13. 71 essentially 
is based on the features of the three last operators in the r.h.s. of (|2.13p we mentioned in Sec. [2j 

Subsequently to Propositions l3.6l and l3.7l we are in a position to isolate a first still _R-dependent 
main contribution from (|3.14p in the tight-binding situation and in the zcro-tcmperaturc regime. 
Under the conditions of (|3.14p . we show that VO < a < 1 there exists a i?- independent constant 
c> s.t. 

\imXp{l3,po,b = 0) = -^XR{£R.F{poiR)),b = 0)+0{e-''^^), (3.16) 

/3too 

with: 

XRi£R,FipoiR)),b = 0):=(^^^ '-j^ de(£fl,F(po(i?))-Ox 

X Tri.(K3){xnJi/p - 0-'[TpA£.)Tp,i{0 - Tp.2(e)]x^^«}. (3.17) 

To achieve that, let Ino be an open interval s.t. /„„ C (A„o,A„o+i) and (A„o + A„o-)-i)/2 G /„(, 
when 7^0 < t; otherwise /„(, C (Ar,0) and Ar/2 G Ina- From Lemma [331 and Proposition 13.11 
it follows that for R sufficiently large /„„ n a{Hp) = and the Fermi energy £fi^p{po{R)) := 

limi^-foo tipHfi, Po{R),b = 0) G /„„. Then p.l6p follows from Propositions 13.61 and 13.71 together. 

Remark 3.8 Ln the case of < t, by virtue of the asymptotic expansion in Lemma \S.l\ iii] along 
with (j3.15p . then one obtains instead of (|3.17l) ; 

lim Xr[P, po{R),b = 0) = -L-XR{£p,F{no),b = 0) + ©(e''^^"), 

0toO \ilR\ 

with XR{£p_p{nQ),b ~ 0) as in (|3.17p but with £p,F(no) defined in p.2p . instead of £r p{pq{R)). 

3.1.3 Isolating the main contribution at zero-temperature - Continuation and end. 

The continuation of the proof of p.2ip consists in removing the i?-dcpcndancc from the leading 
term of p.l7p . arising from the indicator functions xur inside the trace, while keeping a good 
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control on its behavior for large values of R. To achieve that, the crucial ingredient is the expo- 
nential localization of the eigcnfunctions associated with the eigenvalues of Hp in (— oo,0). 



Let us introduce the families {ge.iu, ^ G C}, w ^ 0,1 where Qg^^ : C — > C are defined by: 

3eAO-^0-^, 3edO-=0. (3.18) 

We need the following technical result whose proof lies in Appendix, see Sec. 14.21 

Proposition 3.9 \/d G C, Vw e {0, 1} and V<; e {1, . . . , r} there exist two constants C = C{d) > 
0, c> s.t. 



Vxe 



maxj 



I 

2^ 



d^0e,^(O{(^p-O"'rp,,(e)}(x,x] 



277 



} < C<re- 



c|x| 



J = 1,2, (3.19) 



where Qg,w are the maps defined in p.lSp . and is any positively oriented simple closed contour 
surrounding the smallest eigenvalues of Hp in (— oo,0) while letting outside the rest of the 
spectrum. 

As a result of Proposition [331 getting the trace out the integration w.r.t. ^ and then removing 
the indicator functions in p.l7p . together make appear an additional remainder term which be- 
haves like 0{e^''^) for some i?-independent constant c > 0. In other words, under the conditions 
of (|3.17p there exists a i?-indcpendcnt constant c > s.t.: 



- -Tr, 



X iHp-0-'[Tp.ii0TpA0-TpA0]} + O{e-^''). (3.20) 
Note that due to p.l9p . the leading term in the r.h.s. of p.20p behaves like ©(no) as expected. 

Gathering ((XTC)) . ([XT?)) and (jX^ together, then to complete the proof of PT^ - PT^ . it 
remains to show that the quantity containing the Fermi energy inside the trace of the leading term 
in the r.h.s. of p.20p plays any role (i.e. we can get rid of it without changing the value of the 
trace). This is contained in the below result, whose proof can be found in Appendix, see Sec. 

Proposition 3.10 With the notations of Provosition lKM one has: 

-^Tri.(R3){ / daHp-0'^[Tp,,iOTp,iiO-Tp,2m} =0- 



3.2 Proof of (ii). 

Before starting, let us recall and introduce some notations. Under the assumptions (M')-(-!2K^n)- 
{■s^nd), let {Xi}J^^, r € N* be the set of eigenvalues of Hp in (— oo,0) counting in increasing 
order. Let 7;, / = 1,...,t be positively oriented simple closed contours assumed to be two by 
two disjoint, chosen in such a way that 7; surrounds the eigenvalue A; while letting outside the 
rest of the spectrum of Hp. Let Hp{b), b G M. he the magnetic 'single atom' operator in (|1.18|) . 
From [II Thm. 1.1] (see also [Ml El]), there exists bo > s.t. V|6| < fao, ^J=ili G g{Hp{b)) (the 
resolvent set). Moreover, by virtue of [51 Thm. 6.1], the eigenvalues A/, I = l,...,r are stable 
under the perturbation Hp{b) — Hp for small values of the 6-parameter. Due to the assumption 
{s^nd), then there exists fai > s.t. V|6| < bi, Hp{b) has exactly one and only one eigenvalue 
A; (6) near A/, I = 1, . . . , r which in the first order are given by A; (6) = A; + bei + 0(6), see [35[ 
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Thm. 2.6 in Sec. VIII]. Actually each eigenvalue Xi{- ) can be written in terms of an asymptotic 
power series in b, see e.g. [71 Thm. 1.2]. We denote by Ili{b) the orthogonal projection onto 
the eigenvector corresponding to the eigenvalue A; (6). Gathering all together, then there exists 
< b < min{bo, bi} s.t. V|6| < b each A; (6) lies inside the closed contour ji introduced above. 

We start the proof of ()1.23[) with the following remark. From the Riesz integral formula: 
V|6|<b, \i{b)ni{b) = Hp{b)Ili{b) = ^ f dU{Hp{b) ~ o~\ 1 = 1,. ■.,T. 
Since dimRan(n((fe)) = 1, / = 1, . . . , r by stability of the Afs, then for any ng G {1, . . . , r}: 

no Tio . „ 

V|6| < b, ^ AK6) =Y.^TL-i^m{Hp{b)Mb)] = ^^^L-^i / d^aHpib) - 0^'}- 

(3.21) 

The following result is concerned with the quantity in the r.h.s. of (|3.2ip seen as a function of 
the 6-variable. Recall that for 9 £ C, ge^ui with w = 0, 1 arc the maps defined in p.l8|) . 

Proposition 3.11 There exists a neighborhood X of b = s.t. V0 £ C andyw G {0, 1} the map: 

6^^e,»(6) ^Tri2(R3){/ 9eA0{Hp{b) - (3-22) 

is twice differ entiable on T. Furthermore, its second derivative at b = read as: 
d^^e 



dfe2 



'-{b - 0) := -Ti-mR^){ I diQe.roiOiHp - ir\Tp^mTp,i{i) - Tpaim- (3-23) 



From p.2ip . then the identity (|1.23p straightforwardly follows from the above proposition, 
together with the fact that each of the A;(- )'s is twice difFerentiable in a neighborhood of 6 = 0. 

The rest of this section is devoted to the proof of Proposition 13.111 It is essentially based on 
the so-called gauge invariant magnetic perturbation theory, see [UJ [TSl [TSl lU El IHl 1201 IHl UHl HSl 
[TTl [T^ HS] and references therein for further applications. 

Before starting, let us introduce some notations. Define V(x, y) e the magnetic phase as: 

,/)(x,y):=ie3-(yAx) = -0(y,x), 63:= (0,0,1). (3.24) 

By [m Thm. B.7.2], Vfe G M and £ Q{Hp{b)) the resolvent {Hp{b)-iy^ is an integral operator 
with integral kernel {Hp{b) - O^Ht ) jointly continuous on M'' \ £»,£>:= {(x, y) g R'^ : x = y}. 
Introduce on ^^(M'^) the operators Tpj{b,£,), j = 1,2 via their kernel respectively defined by: 

V(x,y) GR6\i5, rp,i(x,y;6,e) := a(x - y) • (zV, + 6a(x))(i/p(6) - (x, y), 
Tp,2(x,y;6,0 :- ia2(x - y)(i?p(6) - e)-^(x,y). 

From [m Eq. (2.9)] together with [12l Lem. 2.4], Vry > there exists a constant d = 79(77) > 
and a polynomial p{- ) s.t. V6 G M and G C satisfying dist(^, a{Hp{b))) > 77: 

V(x,y)eM6\A |(Hp(6)-e)^'(x,y)|<p(|ei)^^^, '^^^^TT^' ^^-^^^ 

|rp,,-(x,y;&,OI <p(iei)(l + l&l)^ V'" 1^' , J ==1,2. (3.26) 

|x-y| 
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Hence by the Shur-Holmgrcn criterion the operators Tpj{^), j ~ 1,2 arc bounded on L-{ 

\\iHpib)-0-'\\<Pm), \\Tp.,ib,0\\<pi\m + \b\f, (3.27) 
for another polynomial p{- ). For any k e {1, 2}, m e {0, 1} and 6 G M define on M^: 

T??,fc(x,y;6,0:=E(-l)^ E / dzi • • • / dz, (z0(z„ y) - z<^(z„ x))'"x 

X (iJp(&)-e)^^(x,Zi)rp,,,(zi,Z2;&,e)---Tp,,^(z„y;6,e), (3.28) 

where by convention, we set O" = 1. Here i = {ii, . . . , ij} £ {1, 2}-', 1 < j < k and Xj denotes the 
characteristic function defined as: 

.k(i\ — / 1 if *i ^ k 



^ ' [0 otherwise — j — 

Let us note that due to the antisymmetry of </) in p.24p . the terms in the r.h.s. of p.28p containing 
the magnetic phases identically vanish when x = y. Moreover V77 > 0, Vfe S K and € C satisfying 
dist(^, a{Hp{h))) > -q, T™ i,(- , • ; 6, ^) is jointly continuous on R^. This follows by applying j-times 
[H Lem. A.l] together with ([X^ - iP?^ . Furthermore from ((XMI . along with [H 

Lem. A. 2 (ii)], there exists a 6-independent polynomial s.t. for any k £ {1, 2} and m S {0, 1}: 

V(x,y)GR^ |T?,,(x,y;6,OI<p(ICI)(l + l&l)^x| + Zhelltl^ (^.29) 

where, in the case of x 7^ y we have used the rough estimate |0(x, y)| < |y| |x — y| V(x, y) e M®. 
Remark 3.12 In view of (|3^ . we have on M ; 

T°,.i(x,x;6,0 = - / dz(i/p(6)-0-'(x,z)Tp.i(z,x;&,0, 

JR3 



V (x, x; &, = T?,,2(x, x; fe, = - / dz (i/p(fe) - C)"^ (x, z)Tp,2(z, x; 6, 0+ 
fc=i ' ^R-^ 

+ / dzi / dz2(i/p(6)-0"'(x,zi)rp,i(zi,Z2;6,e)rp,i(z2,x;6,C). 

Now let us turn to the proof of Proposition 13. 1 ll It requires two technical intermediary results. 
The first one deals with the regularity of the integral kernel of the resolvent operator (iJp(6) — ^)~^, 
seen as a function of the 6- variable, in a neighborhood of = 0: 

Lemma 3.13 Let K C {g{Hp) n {C £ C ; 5RC < 0}) be a compact subset. Let > s.t. 
V|6| < bK, K C {g{Hp{b)) n {C e C : 5RC < 0}). Then £ K and V(x,y) £R^\D, the map 
b I— >■ {Hp{b) — ^)~^(x, y) is twice differentiable on the interval (— bxj&A')- Furthermore, its first 
two derivatives at bo £ {—bK, bx) read as: 



^(i7p(6)-0-^(x,y) 



b=bn 



= (z0(x, y))^ (iJp(6o) - 0"' (x, y) + sJ2 '^pT/c'Xx, y; 60, « = 1, 2, 



fe=i 



where the functions Sp j,(- , • ; ^Oi = 1, 2, to = 0, 1 are defined in p.28p . 

The proof of the above lemma can be found in Appendix, see Sec. 14.31 The second result 
deals with the 'uniform' exponential localization of the eigenfunctions associated with the A;(6)'s, 
Z = 1, . . . , T for small values of the 6-parameter: 
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Lemma 3.14 Let A;, I G {1,...,t} be a simple eigenvalue of Hp in (— cx3,0). Let b sufficiently 
small s.t. Hp{b) has exactly one and only one eigenvalue Xi(b) near A/. Denote by $/(•;&) the 
associated (normalized) eigenfunction. Then there exists b > and two constants c, C > s.t. 

V|6|<b, VxeM^^ max{|$,(x;6)|,|9^,4>,(x;fo)|} < Ce"'=l''l, fc = l,2,3. (3.30) 

The exponential decay for the $/(-;&)'s is a well-known result, see e.g. Thm. 4.4] and 
also [34j Thm. 1.10], [26l Sec. 7.2]. For b sufficiently small, all the constants can be chosen 
6-independent. 

As a result of Lemma [3. 141 one straightforwardly gets as a corollary of Proposition [ 



Corollary 3.15 Let {\i}J^i be the set of eigenvalues of Hp in (— oo,0). Let b > s.t. V|6| < b; 

(i) . Hp{b) has exactly one and only one eigenvalue Xi{b) located nearby Xi, I ~ 1, . . . ,t. 

(ii) . The (normalized) eigenfunction associated with each Xi(b) obeys the estimate (j3.30|) . 
Then V0 e C, Vw e {0, 1} and V<r G {1, . . . , r} there exist two constants C ^ C{9) > 0, c > s.t. 



V|6| < b, VxeM^, ma^i^ f QeAOiiHpib) - 0''Tp,3ib,0}{^,^) 
^ ^ d^9eA0{{Hp{b) - 0''Tp,,{b, OTpAb, 0}(x, x) 



}<C,{l + \b\)''e-^\-\, J = 1,2, 

(3.31) 

where Qe_w are the maps defined in (|3.18p and F,; the contour as in Proposition \3.9[ 



Let us note that the presence of the factor (1 + \b\)^ in the above upper bound comes from the 
estimate (|3.26|) (we recall that the estimate in (|3.25p is 6- independent). 
We are now ready for: 

Proof of Proposition [Sini Let 6* G C and w e {0, 1}. Let b > s.t. V|6| < b, holds. We 

first prove that: 

V|6| < b, :Fe,Ub) -.^^ f d^Se AO{Hp{b) - 

has an integral kernel jointly continuous on R^. Note that Fg^w{b) is an integral operator since 
{Hp{b) - C)"^ is bounded from X^^rS) i°o(R3) by some polynomial in |^|, see ((3?25l) . Let 
^0 < inf (T(iJp) and large enough s.t. < niiiii<i<r{7i H M}. By the first resolvent equation: 



d^QeAm - iAHp{h) - CrHHpib) - Co)-'- 



I 

2^ 



By virtue of the Cauchy-Goursat theorem, the first term in the above r.h.s. is identically zero. 
The second term has a jointly continuous integral kernel on by virtue of p21 Lem. A.l] along 
with (|3.25p . Denoting by J-g,w{' r \ b) the kernel of J-e.w{b), its diagonal part reads on (— b, b) as: 

VxeM^^ J-e,„(x;6):= J-e,„(x,x;6) = -^f / d^0e,»(O(^fp W - e)"'(x, y)) 

Next from Lemma [3. 131 together with p.29p . one can prove that the map b n- Tg^^ix; b) is twice 
differentiable on (— b, b) Vx G R^. Moreover its first two derivatives at 6o G (— b, b) satisfy: 

VxeK^^ ^lI^[x-bo) := ( de0e,»(e)E^pt(x,x;&o,e), ■'^=1,2. 
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Here we used that 0(x, x) = 0. Finally let < b < b s.t. V|6| < b and for any I £ {l,...,no} 
the (normalized) eigenfunction associated with Xi{b) obeys an estimate of type (|3.30p . From the 
explicit expressions in Rcmark l3. 121 together with the estimate (|3.31l) . then for any compact subset 
K C (— b, b) there exist two constants c > and C = C{no,d, K) > s.t. 

<Ce""l^l, s = l,2. 

The upper bound belonging to i^(R'^), the proposition follows by standard arguments. □ 

3.3 Proof of (iii). 

Let us recall some notations. Under the assumptions (^r)-(-24i)-(-*24d), let {Xi}J^i be the set of 
eigenvalues of Hp in (— oo, 0) counting in increasing order. For any / G {1, . . . , t}, denote by the 
normalized eigenfunction associated with A;, and by Hi = |<I';)($(| the orthogonal projection onto 
the eigenvector We define := 1 — H/. From [SJ Thm. 6.1], we know that there exists bi > 
s.t. V|6| < bi, each A; is stable under the perturbation W{b) Hp{b) — Hp ~ ha- (— iV) + ^^^a^. 
This means that V|6| < bi, Hp{h) has exactly one and only one eigenvalue A; (6) near A; which 
reduces to A; in the limit 6 — >■ 0. For such 6's and any / S {1, . . . , r} denote by ^i{h) the normalized 
eigenfunction associated with \i{h), and by Ili{h) — the orthogonal projection onto 

the eigenvector $;(&). We define n/-(&) := 1 - Xli{b). 

Let us turn to the proof of (iii). Let I £ {1, . . . ,t} and K be a compact neighborhood of A;. 
Let < b < bi s.t. V|6| < b, A/(6) e K and Kn {a{Hp{b))\Xi{b)) = 0. According to the Feshbach 
formula in [22] and under our conditions, V|6| < b A/ (6) is the unique number ^ near A; for which: 

(Az - OUi + UiW{b)Ili - UiW{b){U^{Hp + W{b) - C)Ul}-'Wib)Ui, (3.32) 

is not invertible. Note that 14^(6)11; is bounded on L^(R'^). and its operator norm behaves like 
0(|6|). This follows from the exponential localization of $/ and V$/ in (|4.15p . We mention that 
1^(6)11; remains bounded even with an exponential weight. More precisely, there exists Sq > 
and a C > s.t. 

VO<e<£o, \\W{b)Uic'''-'>\\<C\b\, (•) := + p. (3.33) 

Now let us justify that the operator Il^{Hp{b) - On/- is invertible V|6| < b, G K. To 
achieve that, introduce the Sz-Nagy transformation in [3SJ Sec. L4.6] corresponding to the pair of 
projections n;(6), 11;: 

uib) = (1 - (n;(6) -no^)-^{n,(6)n; + n,^(6)n,^}, 

where Ui{b)Ui -Ul{b)U^ : UiL'^{R^) Ui{b)L'^{R^). Since Ui{b) converges to Ui in norm by the 
asymptotic perturbation theory in |35| . then the above square-root is well-defined by a binomial 
series. U (b) is a unitary operator, and it intertwines both projections: 

Ui{b)^U{b)UiU*{b). 

Note that n/-(6) = Uib)UfU*{b) what implies U*{b)U^{b) = U^U*{b). As a result: 

V|6| < b, e K, U{b)n^{Hp{b) - O^^U*{b) = n^{b)U{b){Hp{b) - OU*mt{b), (3.34) 

i.e. Il^{Hp{b) — is unitary equivalent with the operator in the r.h.s. of p.34p . Next we use 

some results from the asymptotic perturbation theory, see [351 Sec. VIIL2.4]: 

u{b) = 1 + {niWib){Hp - - no - (1 - iii){Hp - i)-^w{b)iii} + ui(5), 

U*{b) = l + {{l-Ui)(Hp-0''W{b)Ui-UiW{b){Hp-0'\l-Ui)}+n2{b), 



Vx G 



sup 



db' 



-(x;6) 
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where the Uj(6)'s, j — 1,2 are operators satisfying b ^Uj{b) in the strong sense. Putting these 
asymptotic expansions into p.34p . then we obtain: 

V|6| < b, G K, um^iHpib) - m^u*ib) = ntmHpib) - mHm + ^it 

where T(^, &) is an operator s.t. |jT(^, — 0{\b\) uniformly in ^ G K. Ergo we eonclude that 
U^{Hp{b) - ^)U^ is invertible V|6| < b and € K. 

Let us get back to the quantity in (|3.32p . By the use of scalar products, the Ai(6)'s has to obey 
the equation: 

v|6| < b, \i{b) = \i + T^(6)$i) - {W{b)<s>i, {n/-[ijp(&) - Xi{b)]ni'}-^w{b)^i). (3.35) 

By iterating the identity in p.35p . one obtains: 

y\b\ < b, Mb) ^Xi + W{b)<i>i) - {W{b)<^i, {UtiHp - Xi)Ut}-^W{b)<i>i) + 0{\b\^). (3.36) 

To control the behavior in b of the remainder term, we used that: 

||e-<->{n^(Fp(6) - A,(&))n^}-^ - {U^iHp X,)Ut}-'c-^^-^\\ = Oi\b\), (3.37) 

with < £ < £0 as in p.33p . p.37p can be proved by using an asymptotic expansion for the 
projection Tli{b) together with the application of the gauge invariant magnetic perturbation theory 
for the kernel of the unperturbed projector as we did for the kernel of the resolvent in Sec. 13.21 
Note that the fourth term involved in the expansion p.36p can be identified with: 

-($i,VK(6)$i)(VK(6)$i,{ni^(ifp-A,)n^}"'VK(6)$0 + 

+ {wmi, {u^{Hp - Xi)ut}-^w{b){n^{Hp - Xi)n^}-'w{b)<i>i). 

Since for b sufficiently small A;(- ) has an asymptotic series expansion, then from p.36p one obtains: 

^{b = 0) = a2$,) - 2(a • (-iV)$i, {U^iHp - A,)n^}-^a • (-zV)$,)- (3-38) 
From (|338l) . the proof of pTM]) directly follows from (fOSl . 



4 Appendix. 

4.1 Proof of Propositions 13.61 and 13.71 

We start by introducing some notations. Recall that under the assumptions (^r)-(=e/m)-('*2^id), 
{^i}J=i with T e N* denotes the set of eigenvalues of Hp in (— oo,0) counting in increasing order. 
For the sake of simplicity, we set A^+i := 0. Let I^, <j € {1,...,t} be an open interval s.t. 
!<; Q {K, K+i) and (A^ + A^+i)/2 G I^. Without loss of generality, we give an explicit form for 
for instance: 

:~ ( _^ll!lj^£±i i±i~) with the convention A^+i := 0. 

For any € {1, . . . , r}, introduce the following decomposition of the contour (recall that 

it is defined as in p.3p but with 6p := inf a{Hp) — 1 = Ai — 1 instead of dn): 

'^r=^Su7Sur.,;3, (4.1) 
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where Lj^ := + 5\^+i)/2A and u!^ {5X^ + 19A5+i)/24; with the convention At-+i := 0. 

Let us start with the proof of Proposition [2111 
Proof of Proposition [nUl Let e {1, . . . , r}. From (|3TT2)) and gj]), V^t G /<,- and Vi? > 0: 

X TvL2^^s){xnAHp - 0-'[Tp,iiOTpAO ~ Tp^^mxn^}- (4.2) 

Since [a;~,a;+] n a{Hp) = 0, then G the closed subset surrounding by 7^^^ is a strict subset 
of the holomorphic domain of the integrand in (|4.2p . Ergo the Cauchy-Goursat theorem yields: 

Now the contours 7^^^ and can be deformed respectively in 7^^^ and r<;^i (set /3 = 1 in their 
definition) due to the location of the interval I^. In the wake of the deformation of contours, then 
under the conditions of Lemma [^TTl from (|2.15l) and p. lip with S = P and by setting = 1, there 
exists a polynomial p{- ) independent of /3 (and R) s.t.: 

|f]«|-i|Tri.(K3){xoJffp - i)-^[Tp,i{OTp.i{0 - rp,2(C)]xn„}| < P(l^l)- (4-3) 
Afterwards it remains to use the Lebesgue's dominated convergence theorem which provides: 

hm i 1^^^ dCf(/?,A^;^)Tri.(M3){xo„(i/p - i)-'[TpAOTp.i{i) - Tp,2{0]xn^} = 0, 

lim i / dai(3,^,;0TrL2^^s^{xnAHp-0-'[TpA0Tp.i{0-TpA0]xn^} = 

f daf^-0'^rLHR^){xnAHp-0''[TpA0Tp.i{0-Tp.2{0]xn^}- 
Here we used the pointwise convergence: for any fixed /i > Ai = inf <7{Hp), 

lim if(/3,pt;^) = (Ai-OX[Ai,M](0- 
The uniform convergence on compact subsets K C /<; is straightforward. □ 
Remark 4.1 FKe emphasize that the deformation of contours "f^^^ andT^_i3 in some P -independent 



contours is crucial in order to make the estimate in (j4.3p f3 -independent. 

Now let us turn to the proof of Proposition 13.71 To do that, introduce some new operators. 
Recall that TZpiO and ^ij(^), R > Ro are respectively defined in (|2.6p and (j2.12p . Introduce 
Vi? > Ro the bounded operators Trj{^) and ^RjiO^ J = 1,2 on L'^{R.^) generated via their 
kernel respectively defined by: 

V(x,y) eM^^^^ rp.i(x,y;e) := a(x - y) • (^Vx)(7^p(0)(x, y), 

rp,2(x,y;e) ia2(x - y)(7^H(e))(x, y), 

=^p,i(x, y; e) := a(x - y) • (jVx)(^fl(0)(x, y), 

^fl,2(x,y;0 := ia2(x - y)(^fl(C))(x, y). 

To prove Proposition 13.71 we need the following four lemmas whose proves can be found in 
Sec. Em Recall that -Ro,i?i > 1 are respectively defined through ([2J|) - p.l4p : 
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Lemma 4.2 Let < a < 1 be fixed. Then for every rj > there exists a constant ~ ^{f]) > 
and a polynomial p{- ) s.t. Vi? > i?o o.'rid G C satisfying dist(^, (T(_ffii) H a{Hp)) > rj: 



1 



1 + ier 



(4.4) 



Lemma 4.3 Let < a < 1 be fixed. Then for every jy > there exists a constant d ~ ^{ri) > 
and a polynomial p{- ) s.t. Vi? > niax{i?o, -Ri} and G C satisfying dist(^, a{HR) D a{Hp)) > ry; 



1 



inn] 



Ti'LHR^){xn^T^RiOrRAOxnn} - T^^-LHm^Axn^^RiO'S^RAOxnJ < 



(4.5) 



Lemma 4.4 Let < a < 1 be fixed. Then for every 77 > there exists a constant ~ ^{rj) > 
and a polynomial p{- ) s.t. Vi? > max{i?o, i?i} and V^ G C satisfying dist(^, a{Hji) D a{Hp)) > r]: 



1 



\^R.\ 



T^rL2(R3){xnRT^Ri07'R,ii0TRs{0xnn}- 

T^LVm^AxnMO^^RAiO-^RAOxnJ < p(|e|)c-^««° 



(4.6) 



Lemma 4.5 Let < a < 1 be fixed. Then for every 77 > there exists a constant d = i9('7) > 
and a polynomial p{- ) s.t. Vi? > i?o and V^ G C satisfying dist(^, CT(ii^) H a{Hp)) > rj: 



1 



TrmM3){xnM0[^RA0^RA0 - ^RAOhnJ- 

Tri2(K3){xoJi?p - 0''[TpAOTpAO - TpAO]xnn} < pdCDc-''^^" . (4.7) 



Now we are ready to prove: 

Proof of Proposition Let i?2 > 1 s.t. Vi? > i?2, M<7{Hr) > Ai - i. A such i?2's exists 
since inf cr(iii^) has to coincide with inf a{Hp) when i? f 00, see Lemma |3^ In view of (|3.8p and 
(|3?T2l) . then V;3 > and V/i G M the quantity Ai?(/3, A) in (IXTSj) can be rewritten Vi? > i?2 as: 



1 1 i 



cj 13 \np] IT j^p 



dCf(/3,/i;0^fl(0, 



(4.8) 



with: 



I^r{0 ■■= Tri.(K3){xn«{(i?fl - 0''[TrAOTrA£.) - TrAC)] + 

- {Hp - 0-'[TpAOTpAO - TpAO]}xnr.}- (4-9) 

Let us estimate (|4.9I) . In view of the definition of the "^j^^-contour, set 77 :— min{l/2, tt/ (2/3)} > 0. 
Let < a < 1 be fixed and define i?3 := max{i?o, i?i, i?2} > 1. From Lemmas 14. 2114. 5[ there exists 
a constant 7? = 7?(7/) > and a polynomial p{ - ) s.t. 



Vi? > i?3, vc G 



AicRm<pm>- 



-^fR" 



l9 



1 + ier 



(4.10) 



Note that the 7?'s and the p{- )'s in (|4.10p are /3-dependent (at least for /3 large enough). From now 
on we limit in (|4.8p the //-domain to the interval i^, <; G {1, . . . , r}, and we suppose that i? > i?3 
and large enough so that [{5X^ + X^+i)/6,{X^ + 5A;;+i)/6] n cr{Hp) = (with the convention 
Ar+i ■— 0). Due to the inclusions C [a;+, w^] C [(5A<; + A.j+i)/6, {X^ + 5A5+i)/6], it is obvious 
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that Vi? > i?3 and large enough, n (j{Hji) = as wcU as 7^"'^,f n (y{Hii) = 0, where 7^"'^ and 

f^^^ are the contours coming from the decomposition of the '^j^'' -contour in (|4.ip . 
Afterwards by mimicking the proof of Proposition I3.6[ then for any compact subset K <Z I^: 

AflX/.) lim Afl(/3,Ai) - f^) / d^ (/i - O/ChIC), 

/3too \c/ |i2/j| TT Jp^^^ 

uniformly in /i € A'. We emphasize that, as in the proof of Proposition 13. 6[ the deformation of 
7^^^ and r<;_^ in some /3-indepcndent contours makes the estimate in (|4.10p /3-indepcndent on these 
contours, see above the definition of the r^'s. Finally it remains to use (|4.4p - (|4.7p along with (|4.10p 
which lead to the existence of two constants c, C > s.t. Vi? > and large enough: 

V^e/,, \^R{^i)\<C{l + \^i\)e-^''\ 

□ 



4.2 Proof of Propositions and I^AOl 

Proof of Proposition 13.91 For any / g {1, ...,?}, let H; be the orthogonal projection onto the 
eigenvector corresponding to the eigenvalue A;. Recall that it is given by a Riesz integral: 

where 7; is any positively oriented closed simple contour surrounding A/ but no other eigenvalue. 
Denote by := ^i- Let us now introduce the following decomposition of the resolvent: 



y^egiHp), {Hp~0'' ^{Hp^O^' + iHp-Ol\ (4.11) 



with: 



{Hp - e), 1 := U^iHp - = E A^^nz, 



(4.12) 



1=1 



{Hp (1 - ^n){Hp = -f / ^(7^(Hp C)-'. (4.13) 



Due to ()4.12p . we have in the kernel sense on 



{Hp - e)^i(x,y) = j^M^)My), (4.14) 



1=1 



where (f>; denotes the normalized eigcnfunction associated with the eigenvalue A/, I = l,...,c. 
Note that under our conditions, the $/'s with ? = 1, . . . ,<r satisfy, see e.g. [lU Eq. (5.8)]: 

VxgM^, max{|$,(x)|,|9^,$/(x)|} < Cc-\/I^I''I < Cc-\/I^I^I, fc = l,2,3, (4.15) 

for some constant C > 0. Moreover denote by Tpj{£_), j = 1,2 and p = H,± the bounded 
operators on L'^{M.^) defined via their kernel as in (|1.19p - (|1.20p but with {Hp — ^)~^ instead of 
{Hp Note that due to (gUS]), the kernel of T^j (0, = 1, 2 still obey the estimate ([5JT|) . 

We are now ready for the actual proof. Under the conditions of the proposition, the first part 
of the proof consists in showing the existence of two constants c > and C = C{9) > s.t. 



Vxe 



< C<re-'=l''l, i = l,2. 
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Let 9 E C By virtue of the decomposition (|4.1ip . introduce the bounded operators on L^(R^): 
-^o'i"''' ^ I d^Qe.UOiHp P = J = h2,w = 0, 1. 

Note that they are integral operators with integral kernel ^g''lj^'^{- , ■ ) jointly continuous on R^. 

At first, due to the location of the F^-contour -^g"^,'^'^ = 0, = 1, 2, w = 0, 1 by virtue of 
the Cauchy-Goursat theorem. Secondly let us look at the diagonal part of the integral kernel of 



^0),^<>J- Pj.qj^ (|4.14p followed by the residue theorem: 



VxeR3, ^^'^'''^i^,^)^y2se,Mi)M^) f dz$,(z)T^.(z,x), j = 1, 2, u; = 0, 1. 

z=i -^R' 

Due to ()4.15p and p. lip , use now that there exists a constant C > s.t. 

Vl<;<^, sup/ dz|¥K^||r^-(z,x)| <C7, J = 1,2. 

Then the /-independent estimate in (|4.15p leads on to: |^g-'^'^'^(x, x)| < Ci^e^^l^'l, j = 1,2, 



w = 0,1 for another C ^ C{6) > and c > 0. Also one has by similar arguments on 



p3. 



^)l — ^'^^ J J ^ 1: 2, ui = 0, 1. Finally, the last terms we have to treat read 
VxeM^, ^,(\)-^'^(x,x) := V $,,(x) / dz¥^a(z-x)•(^V,$^J(z)¥^, 

^,(W(x,x):= ^ $,,(x) / dz¥I>)la2(z-x)$,,(z)¥^, 



as: 



and 



^i5'^^^x,x)=0, .7 = 1,2, 



where we have used the following identity provided by the residue theorem: 

E*^ I 1^ fl9,t«(C) ^ f -2i7r E/i ;2=i when w = 1 

Jr. (A;i-0(A;. -0 " I ' when ^« = 

It remains to use the rough estimate |a(x — y)| < (|x| + |y|) together with the /-independent one 
in (|4J5l) ensuring that Vx e M^, |^g^;'^^'^'^(x, x)| < C?e"'=l''l, j = 1,2 for another C = G{Q) > 
and c > 0. Note that the crucial ingredient involved here is the following estimate: 

Vzy > 0, V^i > 0, re"'^* < Ce"^*, t > 0, (4.16) 

for another constant C = C{v, ji) > 0. 

The second part of the proof consists in showing the existence of two other constants c > 
and C = C{9) > s.t. 



Vx e 



Let 9 € C By virtue of the decomposition (|4.1ip . introduce the bounded operator on L^(M'^): 



-^eT^--^^l d^3eAO{Hp~0^'T^,i{m,iiO, P = N, ±, = 0, 1. 
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Note that it is an integral operator with integral kernel , ■ ) jointly continuous on R^. At 

first ■^Q^^''^ = 0, w = 0, 1 by virtue of the Cauchy-Goursat theorem. Also, by a straightforward 
calculation jy^^"^ = 0, w = 0, 1 since the residue theorem provides us with the identity: 

^ I- (A/1 -0(A/. -0(A/3 -0 



,l2,h = l 



It remains to treat six terms. Let us treat the trickiest one (we make it clear hereafter), that is 
.y^"*^^'^, w = 0,1. From the residue theorem, the diagonal part of its integral kernel reads as: 

VxeM^, ^//'^(x,x) = Vfle,«,(AO / dzi / dz2(i7p-A;)li(x,Zi)x 

X a(zi -Z2) • (iV^i$/)(zi)$;(z2)T^i(z2,x; A,), w = 0,1. 
Let us define VI < Z < and V(x, zi, Z2) G with x 7^ zi 7^ Z2: 

J'/^Hx,zi,Z2) := |(i7p-A0l'(x,zi)(zV,,$0(zi)||zfc||^I(^r^i(z2,x;A0|, k ^ 1,2. 
Start with fc = 1. From (|2.15p . p. lip and ()4.15p there exist two constants C, c > s.t. 



VI < Z < <r, J'/^'(x,zi,z2) < C ^J^^ |z^|e-v^l^ile-%/l^l''^le~'tl^^"^l 



|X-Zi| |Z2-x| 



Using the obvious inequality: "^'"^S'Vi^XI^^I+l^^-xl) < min{f , vTaHIIxI ^ j^s^ ^^j^j^g ^j^j^ 

the uniform estimate obtained from (|4.16p : jzije^Vl'^'^ll^il < |A<;|~^ Vzi G M.^, then there exist 
another constant C > s.t. 

Vl</<c;,VxeR3, [ dzi / dz2 J'/^^(x,Zi,Z2) < Ce-"""^^'V^}I''I. (4.17) 
Here we used that: 



X— zi r ^— ^ Z2— X 



/ e 21- '■ii / 1 c 
sup / dzi -j j- X sup / dz2 -j j- < cste. 

X6K3JR3 |x-Zi| xeR3jR3 |Z2 - x| 

By similar arguments, the upper bound in (|4.17p still holds true when replacing J^/^'' with ij/^''. 
Hence one concludes that Vx G R^, |^;^^'-^(x, x)| < C^e"^!^!, u; 0, 1 for another C C{9) > 
and c > 0. We do not treat the other terms which are simpler. Indeed, they all come from 
operators having the form ^J^^'^ or ^g^^'^ , p =-L,H and w ^ 0,1 which have the peculiar- 
ity that the diagonal part of their integral kernel can always be written on M^, via the residue 
theorem, as ^21=1 with sup^ggs |^i(x)| < c for some constant c > uniformly in I. 

Thus the expected exponential decay in |x| only comes from the /-independent estimate in (j4.15p .D 



Proof of Proposition 13.101 The proof we give requires the notations introduced in Sec. 13.21 
We do not recall them, and refer to the beginning of Sec. 13.21 Let <; G {1, . . . , r}. From the Riesz 
integral formula for orthogonal projections together with the definition p.22p . one has: 

V|6| <b, J2 TrL2(R3){n,(6)} = ^Tr^^^^.^i / d^ {Hp{b) - 0"'} ^- ^1,0(6). 
1=1 •'^Uifi 

Here the b's and the 7i's are the same as the ones appearing in p.2ip . Since the map b 1— ?> ,^i.o{b) is 
twice differentiable in a neighborhood of & = (see Proposition l3. 1 ip . then by using the expression 
for the second derivative at & = given in p.23p (with 9 — 1 and w ~ 0), one gets the identity: 
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But the quantity in the above l.h.s. is identically zero. Indeed by stability of the eigenvalues of 
Hp in (— cx),0), one has diniRann((6) = 1 and thus the sum is a ^-independent quantity (equals 
to ?). □ 



4.3 Proof of Lemma [3A3l 

Let 60 e M. On ^^(rs^ introduee V6 G M and G g{Hp{bo)) the operators ^p(6,6o,0 and 
Tpj(b, &o, C)i J ^ 2 via their kernel respectively defined by: 

V(x, y) e R« \ D, i?p(x, y; 6, 60, := e^^^^^-^^) (ffp(feo) - 0"'(x, y), (4.18) 
fp,,(x,y;6,6o,e) :=c*''^(''''')rp,,(x,y;6o,0, 56:=6-&o, (4.19) 
where <j> stands for the magnetic phase defined in (13.241) . Set also: 

fp(6,6o,C) :=56fp,i(6,&o,e) + ('56)'Tp,2(6,6o,e). (4.20) 

Except for a gauge phase factor, the kernel of Rp{b,bo, £,) and Tp,j{b,bQ,^), j — 1,2 is the same 
as the one of {Hp{bo) — ^)~^ and Tpj{bo,£,) respectively. Therefore V?7 > and G C satisfying 
dist(^, (T{Hp{bo))) > rj, then V& G M Rp{b, bo, ^) and Tp,j{b, bo, are bounded with operator norm 
obeying p.27p (with bg instead of b). Under the same conditions, introduce on L'^{U.^): 

fW(&,6o,0 ■.= -Rpib,bo,Ofp.i{b,bo,0, (4.21) 

T^?) (5, 60, := Rp{b, bo, 0{Tp.i{b, bo, OTp.i{b, bo, - Tp^2{b, bo, 0}, (4.22) 
as well as, with the additional condition G g{Hp{b)) D g{Hp{bo))- 
1 

'I^^\b,bo,0 ■={5bfY.^Sb)'' J2 xl'"'{i)Rpib,bo,OTp.^Ab,bo,OTp,^Ab,bo,0+ 

k=0 iG{l,2}2 

-{Hpib)-^)-\fp{b,bo,Of- (4.23) 

Now we are ready for the actual proof. Let K C [g(Hp) n G C : 5ftC < 0}) be a compact subset. 
From Thm. 1.1], then there exists b^- > s.t. V|5| < bx, K C {g{Hp(b)) n {C G C : SftC < 0}) . 
From now on, let bo G {—bK,bK) be fixed. The starting point of the so-called gauge invariant 
magnetic perturbation theory is the following identity which holds in the bounded operators sense 
on L^{R^), see [HI Proof of Prop. 3.2] and also [H Lem. 3.2]: 

y\b\ < bK,y^eK, {Hp{b) - = Rp{b,bo,0 - {Hp{b) - 0''fp{b,bo,0- (4.24) 

This means that for b sufficiently close to bo, {Hp{b) — ^)~^ can be approximated by Rp(b, bo,£,) 
since the operator norm of the second term in the r.h.s. of (j4.24p behaves like 0(1561). This comes 
from (11201), the definitions (|ITTS)) - (|ITTO)) and the estimates ([X^ - ljX^ yielding: 

y\b\<bK, max{sup|ji?p(6,6o,OII,sup|jfp.j(6,6o,e)||} <C^, J = 1,2, 

for some constant C = C{\bo\,K) > 0. Now by iterating twice ()4.24p and in view of ()4.2ip - (|4.23p : 

2 

V|6| < bK, G K, {Hp{b)-0-' = Rp{b,bo,0 + J2'^^b)''iP{b,bo,0+i^p\b,bo,0- (4.25) 

fe=i 

Afterwards, by rewriting (|4.25p in terms of corresponding integral kernels, one has on K*' \ D: 

2 

{Hp{b) - e)-\x, y) = i?p(x, y; b, bo,0 + 5](<56)'X^'^(x, y; 6, 60, C) + i^i^, T, b, bo,0, (4.26) 

k=l 
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where, for all integer k G {1,2} and for any (x, y) G M^: 

k 

X (i/p(&o) - C)-'(x, zOTp,,, (zi, Z2; 60, • ■ ■ Tp^^, (z, , Y; ^^o, 0, (4-27) 

and Tp ''(•,•; ^0: stands for the kernel of Tp''(6, 6o,0- Let us note that in view of (|4.23p . 
along with and the estimates the kernel T^^(- , • &o,C) behaves like 

0(1(56^) uniformly in ^ G K. Next we remove the 6-dependence in the first two terms of the 
r.h.s. of (|4.26p . To achieve that, wc expand in Taylor power series the exponential phase factor 
appearing in (|4.18p and (|4.27p up to the second order in 5b. Thus V|6| < bif , one gets on \ D: 

Rp{^, y; &, 6o, + jz^5b)^'^f (x, y; b, bo, - ^(<5&)'^^^^^j^(i/p(6o) - O'H^, y)+ 

fc=l fe=0 
2 k 

+ E('5^)' E ~P:r(^'y;^o,C)+4'^(x,y;6,&o,e), (4.28) 

k—l m—1 



where the function Tp , ■ ; ^o, C) is defined in p.28p . and the last term stands for the remainder 

t we have used the explicit expressions in Remark l3.12l to re"w 
term in the r.h.s. of (|4.28p coming from (|4.27p . Note also that by construction, 1p\- ,■ ]b, ba, ^) 



P, 

term. We mention that we have used the explicit expressions in Remark l3.12l to rewrite the second 

(4) 
P 

satisfies the property that its first two derivatives at bo are identically zero. 
Next from the expansion ()4.28p . for b e [— b/^, bx] sufficiently close to bo, it holds on \ D: 

(iJp(6)-e)-i(x,y) - (i/p(6o) -0"'(x,y) = 

(56{*</)(x, y)(i/p(6o) - e)"'(x, y) - ((ffp(&o) - C)-'Tp,i{bo, 0)(x, y)} + o{5b). 
Performing the limit b bo, then the map b i-> {Hp{b) — ^)^^(x, y) is differentiable at bo with 
d 



Q^{Hp{b)-0-\^,y] 



*0(x,y)(i/p(&o) - 0"'(x,y) - mp{bo) - 0^'7^P,i(&o, C))(x, y). 

b=ha 



This result can be extended to the whole of { — bK,bK)- The lemma follows by iterating this 
procedure once again. □ 

4.4 Proof of Lemmas I4.2ti4.5l 

Throughout this section, we denote respectively by || • jj^j and || • the Hilbcrt-Schmidt 
(H-S) norm in 32(i^(K^)) and the trace norm in 3i{L'^{R^)). 

Proof of Lemma SIU Let us denote yR,2{0 |r!priTri2(R3){xo„(i/p - 0"'7'i?.2(0xo«} 
and3;p.i(e) \^R\-^'^^m^^){xnAHR-S.)~^TR^i{OTRAOxn^}- By replacing with 
the r.h.s. of (|2.7p in 3^p.j(^), j = 1,2 then wc have: 

yR.iii) = l^^i^r'Tr,.^(E3){xo„7^p(0rp,l(e)rp,l(C)xa„} + Qfl,i(0, 

where Qr,i{0 and Qp.2(C) consist of seven and three terms respectively. Let 77 > be fixed. 
Let us show that there exists a constant > and a polynomial p{- ) s.t. Vi? > Ro and £ C 
satisfying dist(C, cr(i?p) n cr{Hp)) > 77, |Qp.j(C)| < p(|^|)c"''«-^° j = 1,2. To do that let us take 
some generical terms: 

lR.i(0-=T7^ f dx / dzi / dz2(7^p(0)(x,Zl)a(zl-Z2)•V,,(7^p(0)(zl,Z2)x 

X a(z2 - x) • V,,{(ifp - 0~'Wp(C)}(z2,x), 



30 



qB.,2iO ■■= -T7^ [ dx / dz(7ei?(e))(x,z)ia2(z-x){(if^-e)-^Wi?(e)}(z,x). 

Now we need the following estimates. From (j2.15p . p.l6p and p.20p . there exists a constant > 
and a polynomial p{- ) s.t. Vi? > Rq and G C satisfying dist(^, <j{Hfi) n a{Hp)) > i]: 



V(x,y) e]R^ 
V(x,y) eM*^\A 



/ dz(iJ;^-0~'(x,z)(WH(0)(z,y) 
/ dzVx(i/fl-e)"'(x,z)(W«(C))(z,y 

JE3 



<p(iei)e"''«^°e-^«l"-y|, (4.29) 



<p(l^l)e-'^^^° '[J'l^y|' - (4.30) 



Then from (PT7|) . and g^Hl-dlSni) together with [H Lem. A.2 (ii)], max{|to,i(^)|, |(7fl,2(OI} < 

p(|^|)e~''«^ for another constant > and polynomial p(- ) both i?- independent. 

The others terms coming from Q_rj (C), j = 1,2 can be treated by using similar arguments. □ 

Proof of Lemma HH Let us consider ^i?,2(0 := \^R\~^^T^L^R^){xnnT^R{OTRMOxnn}- By 
replacing TZr{£,) with the r.h.s. of (|2.8p in ?Vr^2{0: then we have: 

?^rAO = I^^Hr'Tri.(R3){xa„^i?(0-^i?,2(aXOH} + ^ii,2(0, 

where ^i?,2(C) consists of three terms. Let 77 > be fixed. Let us show that there exists a constant 
■d > and a polynomial p{-) s.t. Vi? > max{i?o,i?i} and € C satisfying dist(^, cr(i7_R) n 
(^{Hp)) > T], \^R,2{0\ < P(l^l)e"'''=^°- To do that let us take a generical term of ^r,2{0- 

qfl,2(e):--^/ dx/ dz(^^,(C))(x,z)la2(z~x)(#i^X0)(z,x). 

Let us introduce the operators ^i?(C), t|^^(0 and T^\i) = Ts,2{£.): with S := i? or P and 
R 1^ Ro generated via their kernel respectively defined by: 

V(x,y) eM^ irfi(x,y;C) := ia2(x - y)(5rH(0)(x, y), 
V(x,y)GM6\i?, rW(x,y;e) :-a(x-y)(i7H-e)-'(x,y), (4.31) 
rf (x,y;0 - TH,2(x,y;0 :- ia2(x- y)(i/H - e)"'(x,y). (4.32) 

Due to the estimates p.l9p and p.lSp . all these operators are bounded on L^(R'^). Moreover 
from (|2.17p and p.l9p . the operator ^r{£^)S'r{£,) has a jointly continuous kernel on K^, see [l2l 
Lem. A.l]. Let us now prove the existence of a constant d > and a polynomial p{-) s.t. 
Vi? > max{i?o, Ri} and e C satisfying dist(^, <7{Hr) n a{Hp)) > 77: 

WxnAiO^RiOxnnK <p(ICI)e-''^^°. (4.33) 
To do that, use that i^(^) can be rewritten as: 

STRiO - Ir{T^rH.OVr{Hp - + {Hr - 0-'VrTP{0 + T^^\OVRT'it\S,m - 9r), 

where we used that a'^(x — y) = {a(x — z) + a(z — y)}^ for any x, y, z G M"^. Now it remains to 
use these estimates on H-S norms which hold for any R > max{i?Oj Ri}'- 

iTia^{\\xnnMiHE-0-'h,A\xnjR{Hs-0''hA<PmR~\ (4-34) 
m^x{\\VR{Hs - 0"'(1 - 9R)xn„h,, || VhT|^'^(0(1 " 9R)xnj3,} < p(|e|)i?^e-''««° , (4.35) 

with S = i? or P and for another i?-independent polynomial p{-). (|4.34p - (|4.35p are obtained 
from the definitions (|L3T|) - (|L32)) together with ((2?T5l) and Lemma [2^ (iii). Thus \qR,2iO\ < 
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p(|^|)e ^i^" , The two others terms of ^r^2{£,) can be treated by using similar arguments. □ 

Proof of Lemma [1:1 Let us denote: &^rj{0 \^R\~'TrL2^K3){xn^nR{0TRA0TR,i{0xnJ. 
By replaeing TIr{£,) with the r.h.s. of (|2.8p in then we have: 

= \nRr'TTL.(^s){xnMO^RAO^RAOxn^} + ^rAO, 

where J2r^i{£,) eonsists of seven terms. Let ry > be fixed. Let us show that there exists a constant 
1? > and a polynomial p{-) s.t. Vi? > max{i?o,-Ri} and G C satisfying dist(^, cr(_ffi^) n 
CT(-ffp)) > V, \^R,i{0\ < Pd^De"''^^"- To do that, take a generieal term from ^R,iiO- 

^rAO-^tt^ [ dx / dzi / dz2(^fl(0)(x,zi)a(zi -Z2) • Vz,(^,?(0)(zi,Z2)x 

xa(z2-x).V,,(#k(e))(z2,x). 

Let us note that from (^3)) and (PTTU|) : 

V^fl(0 = [(V5fl)(i/P - 0"' + .9i?V(i/p - 0"'].9P + [(V|flJ(i7p - e)-' + ,^flV(i7p - - ffp), 

"^^RiO = [(V|p)(i/p - C)''^i?.(^P - 0"' + .^pV(i/p - 0~'Vr{Hp - - gp). 

This means that by using the definitions (|4.3ip . the quantity qfl.i(^) ean be rewritten as: 

qi?,i(0 = qka(0 + qp,\(e), where: 



qk^.\(e) 





dzi / 


Inn 


S3 Jh 



(V|,j)(zi)rW(zi,Z2;e)(l-gp)(z2)}a(z2-x).V,,(#k(e))(z2,x), 



^rAO - ^^^^LHm{xnn^R{0[9RTpA09R + .^prp,i(e)(l - .9p)]x 



X 



{{ylR)[T'i\S,)VR{Hp + {Hr 0-'VrT^p\0 



Ir[TrAOVr{Hp - + y{HR - 0"'VHT(,i)(e)]}(l - gp)xo«}. (4.36) 



On the one hand due to (|2.3p and p.Sp , by mimicking the proof of p.2ip , there exists a constant 
1? > and a polynomial p{- ) s.t. Vi? > Rq and G C satisfying dist(^, (j{Hr) n a{Hp)) > 77, 

max{|(Vgp)(x)4i)(x,y;0gpXy)M(V|p)(x)4i^(x,y;0(l-.9pJ(y)|}<Mie^ 
On the other hand, from (|2.15p . (|2.16p and our assumption on u: 



Vi?>i?o,V(x,y)GR<'\i?, |Vx(#'p(0)(x,y)| <P(|^|)- 



|x-y| 



for another i?-independent d > and p(-). Hence (|2.17p along with the two above estimates 

lead Vi? > i?o to |qp_\(C)l < p(ICI)e"'''=-^°, see e.g. [H Lem. A.2]. Moreover, by mimicking the 
proof of (|4.33p . the trace norm of the operator inside the braces in (|4.36p is bounded above by 
polynomial x R^c~^^^ Vi? > max{i?o,i?i} due to the H-S norms (|4.34p - (|4.35p and the operator 
norms in ([2:221) - ([2:231) . It follows that Vi? > max{i?o, i?i}, |qp.\(OI < P(ICI)e"''«^° for another 
i?- independent > and polynomial p{-). Hence |qp,i(C)| < Pd^De"''^^" ■ The others terms 
coming from =Sp_i(^) can be treated by using similar arguments. □ 
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Proof of Lemmaimi Let us denote: 2)i?,2(C) := \^r\ ^Tri2(R3){xnfl^jj(f)^,2(C)xo„ } and 
2)i?,i(?) := \^R\~^T^^LHR^){xnn^RiO.^RAO-^RAOxnJ- By replacing ^^(0 with the r.h.s. of 
(Pl^ in 2)/?j (C), J = 1, 2 then we have:: 

^RAO^\^Rr'^^L^M^){xnAHp-0''TpAOTp.i{Oxnn}+QRAO, 
^rAO = \nn\-'TrL-2^^s){xn,,iHp-0''TpA0xnJ + QRA0, 

where 13_r,i(C) and Oi?,,2(0 consist of seven and three terms respectively. Let > be fixed. 
Let us show that there exists a constant > and a polynomial p{- ) s.t. Vi? > i?o and S C 
satisfying dist(^, cr(i?ij) n cr{Hp)) > ?y, < p(ICI)c"''«-^° j = 1,2. To do that let us take 

some generical terms: 

<^rAO--^t^J dx/ dz(ifp-0-'(x,z)ia2(z-x)(2n;^(C))(z,x), 



^rAO 




X a(zi -Z2) • Vzi(ilp -0"^(zi>Z2)a(z2 -x) • Vz,(2Ufl(^))(z2,x). 

From (HHni), (PTTB)) . (P?^ together with [HI Lem. A. 2 (ii)], then we straightforwardly get the 
existence of a constant > and a polynomial p(- ) s.t. Vi? > i?o and S C satisfying 
dist(^, cr(7Jp) n<7{Hp)) > ?7, |qpj(C)l < p(l^l)c"''«^° j = 1,2. The other terms coming from 
Qpj (^), j — 1,2 can be treated by using similar arguments. □ 
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